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1.  INTRODUCTION  AND  SUMMARY 


This  report  presents  major  accomplishments  and  signif¬ 
icant  findings  of  atmospheric  research  performed  under  AFGL 
Contract  F19628-81-C-0016  for  the  period  December  1980  - 
June  1982.  Most  of  the  findings  presented  here  are  based 
on  the  AER  one-  and  two-dimensional  models,  although  selec¬ 
tive  findings  from  other  modeling  groups  are  also  included 
for  comparison  purposes.  Part  of  the  modeling  effort  is 
also  funded  by  the  CMA  modeling  program. 

This  report  is  organized  into  five  sections  together 
with  two  supporting  appendices.  Section  2  gives  a  brief 
overview  of  stratospheric  chemistry  with  emphasis  on  the 
impact  of  recent  rate  data  revisions  on  the  quantitative 
understanding  of  various  stratospheric  chemical  processes. 

The  new  rates  for  the  reactions  of  OH  with  HN0o  and  HNO. 

3  4 

introduce  additional  coupling  between  HO  and  NO  species 
with  the  effect  of  reducing  the  OH  and  H02  concentrations 
in  the  lower  stratosphere.  Calculations  are  shown  to  illus 
trate  the  impact  of  the  new  rates  on  the  vertical  distribu¬ 
tions  of  those  trace  gases  that  participate  either  directly 
or  indirectly  in  various  stratospheric  03  removal  processes 
Section  3  presents  the  results  on  the  potential  ozone 
perturbations  due  to  various  human  activities.  It  is  found 
that  mainly  because  of  the  rate  data  revision,  the  calcu¬ 
lated  steady  state  0^  depletion  due  to  continued  release  of 
,CFCs  at  1977  production  level  is  6.1  percent,  compared  to 


the  previous  depletion  calculation  of  16.4  percent  based  on 
the  NASA  (1979)  rates.  Sensitivity  studies  of  03  perturbed 
by  increasing  N20  and  C02  have  also  been  performed.  It  is 
pointed  out  that  03  perturbations  due  to  various  human  activ¬ 
ities  may  interact  in  a  non-linear  fashion.  Increase  in 
atmospheric  C02  for  instance#  may  partially  offset  the  CFC's- 
induced  03  perturbation,  particularly  on  an  intermediate  time 
scale  (year  1980-2010). 

Section  4  presents  the  results  from  our  2-D  model  cal¬ 
culation.  We  began  with  a  comprehensive  comparison  of  the 
calculated  trace  gas  distribution  with  observations.  In 
general#  good  agreement  is  found  between  calculated  and 
observed  distribution  of  long-lived  species,  except  for 
CFC-11  and  12.  The  calculated  latitudinal  and  seasonal 
column  03  agrees  well  with  observations#  although  the  cal¬ 
culated  latitudinal  and  seasonal  gradients  in  column  03  are 
less  pronounced  than  those  observed. 

Much  better  agreement  between  calculated  and  observed 
distributions  of  radical  species  in  the  lower  stratosphere 
(15-30  km)  is  achieved  when  the  new  kinetic  rate  data  set 
is  used.  Agreement  between  calculated  and  observed  distri¬ 
butions  of  Cl#  CIO,  NO  and  HN03  in  the  35-40  km  region#  how¬ 
ever,  is  rather  poor.  The  situation  is  worth  noting,  since 
the  discrepancies  occur  in  the  region  where  the  largest 
impact  of  CFCs  on  03  is  predicted;  where  the  uncertainty  in 
transport  is  thought  to  have  minimal  effect  on  model  calcu¬ 
lation,  and  where  the  chemistry  appears  to  be  relatively 


simple  because  of  the  reduced  role  of  larger  molecules  {e.g., 
C1N03,  HN04,  HN03)  in  stratospheric  coupling  processes. 

We  also  calculated  the  atmospheric  lifetimes  for  ^0, 
CFC-11  and  CFC-12  using  the  2-D  model.  The  behavior  of 
these  gases  is  typical  of  the  upward  diffusing  gases  which 
enter  the  stratosphere  by  upward  diffusion  from  the  tropo¬ 
sphere  and  are  removed  in  the  stratosphere  by  photolytic 
reactions.  Although  direct  comparison  of  the  2-D  and  1-D 
results  is  difficult,  our  analysis  points  to  the  problem 
associated  with  the  interpretation  of  the  1-D  model  calcu¬ 
lated  lifetimes.  Because  of  the  inability  of  the  1-D  models 
to  incorporate  latitudinal  covariance  effect,  1-D  models 
tend  to  overestimate  the  lifetimes  of  upward  diffusing  trace 
gases. 

In  Section  4  we  also  summarized  the  areas  where  signif¬ 
icant  improvement  has  been  made  to  the  2-D  model.  These 
include  the  development  of  the  diurnal  code  and  a  semi- 
implicit  scheme  for  solving  the  trace  gas  equation.  Other 
ways  for  improving  the  present  model  will  also  be  discussed. 

Recently,  interests  on  meteoric  metals  in  the  atmosphere 
were  renewed  because  of  the  possibility  that  trace  metals 
may  react  with  the  chlorine  species  in  the  atmosphere.  In 
Section  5  we  report  on  our  modeling  study  of  the  subject. 

We  conclude  that  metal  species  have  the  potential  to  couple 
to  other  chemical  cycles  and  change  radical  species  concen¬ 
trations.  However,  quantification  of  its  effect  remains 
difficult  because  of  the  lack  of  kinetic  data.  Nevertheless, 
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further  modeling  wor’^  should  be  helpful  in  identifying  the 
role  of  the  reactions  and  suggest  priority  for  their  measure¬ 
ments. 
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2.  STRATOSPHERIC  CHEMISTRY:  GENERAL  REVIEW 

This  section  is  intended  to  review  the  recent  develop¬ 
ments  in  stratospheric  chemistry.  We  shall  attempt  in  what 
follows  to  provide  a  simple  account  of  stratospheric  chem¬ 
istry  and  focus  on  several  recent  revisions  in  kinetic  rate 
data  that  have  changed  our  quantitative  understanding  of 
several  key  chemical  cycles  responsible  for  ozone  removal. 

Ozone  is  formed  by  photolysis  of  Oj,  followed  by  the 
three-body  reaction 

o  +O2+M  ■+•  Oj  +  M  (2.1) 

Loss  of  O^  is  affected  by  processes  such  as: 

0  +  03  ♦  202  (I) 

and 

03  +  03  ♦  302  .  (II) 

Both  processes  (I)  and  (II)  may  be  catalyzed  by  a  complex 
suite  of  reactions  involving  a  variety  of  trace  gases.  Key 
examples  for  process  (I)  are: 


and 


NO  +  03  ♦  N02  +  02 

N02  +0  ♦  NO  +  02 


(1.1) 
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Cl  +  03  -*■  CIO  +  02 

CIO  +  0  -*■  Cl  +  02 


(1.2) 


The  efficiency  of  catalytic  cycles  (1.1)  and  (1.2) 
depends  on  the  availability  of  0  atoms  as  well  as  the  rad¬ 
icals  NO  and  Cl.  Since  the  concentration  of  0  decreases 
rapidly  with  decreasing  altitude,  the  importance  of  pro¬ 
cesses  such  as  (1.1)  and  (1.2)  as  sinks  for  03  is  often 
restricted  to  altitudes  above  25  km.  The  availability  of 
stratospheric  NO  and  Cl  is,  in  part,  limited  by  the  supply 
of  source  molecules  (e.g.,  N20,  CFCs)  from  the  troposphere, 
and,  in  part,  regulated  by  H0X  chemistry,  an  area  that  has 
been  affected  most  by  recent  changes  in  rate  data. 

In  contrast  to  process  (I) ,  recombination  of  03  with  03 
(process  II)  does  not  depend  on  0  atoms,  but  either  involves 
slower  reactions  which  are  difficult  to  study  at  the  labora¬ 
tory  or  a  longer  sequence  of  reactions.  Key  examples  of  (II) 
are: 


and 


OH  +  03  ■*  H02  +  02 

H°2  +  03  -*•  OH  +  202 

Cl  +  o3  -*•  CIO  +  o2 
CIO  +  H02  -*■  H0C1  +  0 

H0C1  +  hv  -*■  OH  +  Cl 
OH  +  03  -*•  H02  +  02 


(II. 1) 


(II. 2) 
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According  to  our  present  knowledge,  the  N0X  catalytic 
<:y<;  i  *•  (1.1)  1h  tho  major  wink,  tor  stratospheric  in  the 
region  2i>-36  kin.  Removal  by  (1.2)  constitutes  less  than 
10%  of  (1.1)  below  30  km  but  approaches  20%  of  (1.1)  at 
around  40  km.  One  major  concern  is  the  possibility  that 
this  sink  (1.2)  for  O3  might  grow  in  time  as  the  CFCs  con-' 
tinue  to  increase  in  the  atmosphere. 

Ozone  removal  by  (II. 1)  and  (II. 2)  is  restricted  in  the 
lower  stratosphere  and  depends  very  much  on  the  availability 
of  HOx  radicals.  These  radicals  are  mainly  produced  by  the 
reaction  of  0(^0)  with  H20, 


O^D)  +  H20  «*•  20H  ,  (2.2) 

X 

and  are  removed  by  the  following  reactions: 


OH  +  HNQ3 

-* 

H20 

+ 

NO  3 

(2.3) 

OH  +  HN04 

•* 

h2o 

products 

(2.4) 

OH  +  H202 

h20 

+ 

HQ2 

(2.5) 

* 

OH  +  H02 

’  •+ 

h20 

+ 

°2  ‘ 

(2.6) 

Prior  to  1979,  reaction  (2.6)  was  thought  to  be  the 
major  sink  for  stratospheric  H0X.  Because  of  the  complete 
lack  of  data  on  OH  and  H02  below  30  km,  there  is  no  direct 
way  to  test  model  predictions  of  H0„  in  this  region.  In- 
stead,  Sze  (1978,  1980)  and  McConnell  and  Evans  (1978)  have 
chosen  an  indirect  approach  to  validate  HOx  models.  They 
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noted  that  certain  model  features  (such  as  gradients  of  CIO 
and  C2Hg,  HF/HC1,  HNO-j/NO^  would  agree  well  with  observa¬ 
tions  if  it  was  assumed  that  the  models  (prior  to  1979)  had 
overestimated  the  OH  abundance  below  30  km.  New  laboratory 
measurements  lend  support  to  the  "low  OH"  hypothesis.  Spec¬ 
ifically,  the  new  rates  for  (2.3)  and  (2.4)  are  found  to  be 
much  faster  than  formerly  believed  (see  Table  2-1) .  The 
dominant  sinks  for  HOx  below  30  km  are  now  effected  by  (2.3) 
and  (2.4)  and  the  importance  of  (2.6)  is  restricted  alti¬ 
tudes  above  30  km,  as  indicated  in  Figure  2-1. 

Mainly  as  a  result  of  new  laboratory  studies  of  ^) , 
(2.4),  and  of  the  cross-section  of  H02N02,  the  calculc_ed 
altitude  distributions  of  many  key  radical  species  have 
undergone  significant  changes  of  the  past  18  months.  Some 
of  the  examples  based  on  the  AER  one-dimensional  model  are 
shown  in  Figures  2-2  through  2-7.  Note  that  the  OH  concen¬ 
tration  is  about  a  factor  of  2  to  3  smaller  than  that  pre¬ 
dicted  by  1979  models  for  reasons  discussed  above.  A  some¬ 
what  larger  reduction  is  calculated  for  H02  (Figure  2-3) . 

Figures  2-4  and  2-5  show  the  calculated  changes  in 

chlorine  species  (CIO,  HC1,  HOC1,  CINO^)  in  the  stratosphere. 
Many  of  the  changes  are  related  to  the  reduced  stratospheric 
abundances  of  OH  and  H02  predicted  by  the  present  model.  The 
concentration  of  CIO  with  1981  rate  data  is  lower  than  that  of 
the  1979  model  by  factors  of  about  ten  and  two  at  20  and  30 
km  respectively.  The  calculated  changes  for  HC1  (Figure  2-4) 


TABLE  2-1 


Rate  Constants  for  Key  H0X  Reactions  Given  in 
NASA  (1979)  vs.  NASA/WMO  (1981) 


Reaction 


OH  +  HN03  -*•  product 
OH  +  HO2NO2  -*■  product 
OH  +  H02  *  H20  +  02 
OH  +  H202  -*•  H20  +  H02 


Rate 

Constant 

NASA  (1979) 

NASA/WMO 

(1981) 

8.5 (-14) 

1.5  (-14)  exp(^; 

5.0  (-13) 

4  (-12) 

3. 0 (-11) 

8.0  (-11) 

1.0  (-11)  expr-^J  2 . 7 (-?  2)  exp(”^) 


ALTITUDE  (km) 


« 


[t 


Figure  2-1 

Stratospheric  budget  for  odd  hydrogen  radicals  showing 
the  production  and  removal  rates  due  to  various  chem¬ 
ical  reactions.  The  calculation  was  performed  using 
the  reaction  rate  constants  from  JPL  81-3  (1981) . 


13 


he 


»  ■■  *  1  1  *  * 

KT 

OH  DENSITY  (cm’5) 


■  ■  ■  ■  i  » i 

10 


P) 


Figure  2-2 

Comparison  of  1-0  model  calculated  OH  concentrations 
using  reaction  rate  constants  from  NASA  (1979)  and 
NASA/WMO  (1981) 
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» 


ALTITUDE  (km) 


Comparisons  of  1-D  model  calculated  CIO  and  HC1  concentrations  using 
reaction  rate  constants  from  NASA  (1979)  and  NASA/WMO  (1981). 


are  less  dramatic,  since  HC1  is  a  major  component  of  C1X, 
which  is  mainly  controlled  by  atmospheric  transport  and  the 
supply  of  Cl-source  molecules  such  as  CFCs,  CH^CCl^,  CCl^ 

and  CH3C1  from  the  troposphere. 

Major  changes  also  are  found  in  the  calculated  C1N03 
and  H0C1  concentrations.  The  much  lower  value  calculated 
for  C1N03  (Figure  2-5)  is  mainly  attributed  to  the  use  of 
a  slower  formation  rate  for  C1N03  as  recommended  recently 
by  NASA/WMO  (1981) .  It  should  be  noted  that  there  are  still 
considerable  uncertainties  regarding  the  C1N03  chemistry. 
Further  laboratory  studies  together  with  atmospheric  measure¬ 
ments  are  needed  to  resolve  the  issue.  The  largest  calcu¬ 
lated  change  amon  all  species  occurs  in  H0C1,  a  molecule 
formed  by  the  reaction  of  H02  with  CIO.  The  more  than  a 
factor  of  10  reduction  in  HOC1  concentration  (Figure  2-5) 
around  20  km  reflects  lower  concentration  of  both  CIO  and 
H02  in  the  stratosphere.  As  a  result,  H0C1  should  play  a 

much  lesser  role  in  03  removal  through  (II. 2). 

Data  on  HOCl  are  currently  non-existent.  Since  HOC1 

is  sensitive  to  relatively  small  changes  in  rate  constants, 
a  detection  of  this  molecule  should  provide  an  important 
check  to  current  models. 

Figure  2-6  shows  the  calculated  NO  and  N02  profiles 
calculated  with  1979  and  1981  rate  data.  The  model  with 
the  recent  rate  data  predicts  significnatly  higher  NO  and 
N02  concentrations  in  the  stratosphere  (z  <  30  km) ,  where 
OH  is  a  major  sink  for  N02  through  the  reaction  of  OH  +  N02. 


17 


ALTITUDE  (km) 
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Figure  2-5 

Comparison  of  1-D  model  calculated  C1N03  and  HOC1  using 
reaction  rates  from  NASA  (1979)  and  NASA/WMO  (1981) . 
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The  higher  concentrations  of  NO  and  NO2  would  imply  a  more 
efficient  N0X  cycle  (1.1)  in  ozone  removal. 

Figure  2-7  shows  the  vertical  profiles  of  calcu¬ 

lated  with  1979  and  1981  rate  data.  The  observed  vertical 
gradient  of  CjHg  was  identified  as  a  sensitive  indicator 
for  stratospheric  OH  and  Cl  (Rudolph  et  al.,  198;  Sze  et 
al . ,  1981).  Our  models  with  1979  rate  data  predict  a  much 
steeper  gradient  than  observations.  The  1981  model  which 
calculates  lower  stratospheric  OH  and  Cl  predicts  substan¬ 
tially  higher  C2Hg  concentration  around  20-25  km,  although 
the  new  model  profile  above  25  km  still  lies  below  observa¬ 
tions  by  factors  of  3  to  6.  Further  improvement  between 
observation  and  model  is  possible  by  further  reducing  the 
stratospheric  OH  and  Cl  concentrations  predicted  by  the 
current  models.  Alternately,  a  faster  vertical  eddy  dif¬ 
fusion  coefficient  may  also  improve  the  agreement  although 
this  approach  would  further  degrade  the  already  poor  agree¬ 
ment  between  observed  and  calculated  CFC-11,  as  discussed 
in  Section  4. 

The  calculated  changes  in  the  stratospheric  distribu¬ 
tions  of  the  above  species  have  several  important  implica¬ 
tions  for  atmospheric  models.  First,  in  regards  to  the 
process  of  comparing  observations  with  models,  a  change  in 
model  results  can  either  improve  the  agreement  or  introduce 
additional  discrepancies.  On  balance,  the  profiles  predicted 
by  models  with  1981  chemistry  agree  with  observation  better 
than  earlier  calculations  (NASA,  1979) ,  though  difficulties 
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>n  or  1-D  model  calculated  02%  concentrations  using  different 
rate  constants  taken  from  NAS?.  (1979)and  NASA/WMO  (1981). 


still  remain,  particularly  in  the  upper  stratosphere  (MO  km) 
A  detailed  comparison  of  models  with  observations  is  given  in 
Section  4. 

The  second  implication  concerns  the  chemistries  of  cer¬ 
tain  large  molecules  (C1N03,  H0C1,  HNO^  and  HN04)  formed  by 
various  coupling  processes.  The  use  of  new  rates  reduces 
the  roles  of  C1N03  and  H0C1  in  straot spheric  chemistry  but 
enhances  those  of  HNO^  and  HNO^  in  the  removal  of  HOx. 

The  third  implication  concerns  models'  sensitivity  to 
perturbation  studies.  This  issue  will  be  addressed  in  de¬ 
tail  in  Section  3.  Briefly  stated,  the  recent  revisions  in 
rate  data  reduce  the  potential  impact  of  CFCs  but  increase 
that  of  N20  on  stratospheric  ozone. 
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3.  OZONE  PERTURBATIONS  STUDIES 


3.1  Single  Perturbation  Studies 

There  are  various  man' s  activities  that  may  potentially 
alter  the  global  distributions  of  atmospheric  ozone.  Among 
the  more  important  ones  that  have  been  identified  are: 

1)  Injection  of  N0X  from  hypothetical  high  flying 
supersonic  aircraft. 

2)  Release  of  CFCs  due  to  their  widespread  use  as  a 
refrigerant  and  plastic  foam  manufacturing. 

3)  Release  of  NjO  due  to  the  application  of  nitrogen 

fertilizers. 

4)  Increase  of  C02  due  to  the  combustion  of  fossil  fuel 

5)  Injection  of  N0X  from  subsonic  aircraft  operating 
in  the  upper  troposphere. 

6)  Increase  of  CH^  due  to  agricultural  practice  and 
tropospheric  OH  reduction. 

In  addition  to  the  above  possible  anthropogenic  influ¬ 
ences,  ozone  also  may  be  perturbed  by  a  variety  of  natural 
phenomena  such  as  volcanic  eruptions,  solar  proton  events, 
variability  of  the  solar  flux  and  the  water  vapor  in  the 
stratosphere.  The  quantitative  effects  of  these  natural 
events  on  0^  is  very  uncertain  and  their  assessment  is 
beyond  the  scope  of  the  present  report. 

Table  3-1  presents  the  calculated  perturbations  to 


ozone  due  to  various  man's  activities.  Also  included  in 
Table  3-1  are  the  corresponding  model  results  based  on  the 


TABLE  3-1 


id 


Calculated  Change  in  Column  03 
due  to  Individual  Perturbation 


Change 

in  Column  0^  (%) 

L 

Scenario 

1979* 

1981** 

1. 

NOx  injection  of  2  x 10® 

+  4.0 

-  6.2 

at  20  km 

w 

2. 

CFC's  release  at  1977 

production  rates 

-16.4 

-  6.1 

3. 

Doubling  N20  flux 

+  1.6 

-  9.5 

X 

4. 

Doubling  CO  +  narrow  band  model 
y  2  wide  band  model 

- 

+  2.5 
+  3.9 

*  Based  on  NASA  (1979)  rate  data. 

**  Based  on  NASA/WMO  (1981)  rate  data. 

t  The  narrow  band  and  wide  band  models  refer  to 
different  treatments  of  radiation  transport 
calculations  in  two  separate  model  simulations 
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earlier  rate  data  recommendations  (NASA,  1979).  For  the 
CFC's  constant  release  scenario,  we  summarize  the  impact  of 
individual  revision  in  rate  data  on  ozone  depletion  calcula¬ 
tion  in  Table  3-2. 

As  discussed  in  Section  2,  the  recent  revisions  in  rate 
constants  for  several  H0X  reactions  have  changed  the  quanti¬ 
tative  understanding  of  several  key  chemical  cycles  that 
control  the  removal  rates  of  stratospheric  ozone.  In  brief, 
the  reduced  stratospheric  OH  level  calculated  by  current 
models,  implies  a  less  efficient  H0x  and  C1X  cycle  but  a 
more  efficient  N0X  cycle  for  ozone  removal,  particularly  in 
the  lower  stratosphere  (20-30  km) . 

The  calculated  steady  state  column  ozone  reduction  is 
now  6.1%  compared  to  the  corresponding  figure  16.4%  based 
on  the  earlier  NASA  (1979)  rate  data.  Much  of  the  drastic 
change  in  the  depletion  calculation  occurs  below  32  km  as 
indicated  in  Figure  3-1.  Note  that  a  small  increase  in  03 
around  20  km  is  predicted  by  the  1981  model.  It  should, 
however,  be  noted  that  the  new  rate  data  have  relatively 
little  impact  on  the  depletion  calculation  above  32  km. 

The  largest  percentage  change  in  local  ozone  under  steady 
state  condition  is  predicted  to  occur  at  around  40  km  with 
a  magnitude  of  about  42%.  The  calculated  changes  in  column 
03  on  an  intermediate  time  scale  (year  1970-2000)  are  given 
in  Figure  3-2.  The  predicted  present  day  (1980)  03  deple¬ 
tion  is  0.65%,  compared  to  the  corresponding  value  of  1.7% 
predicted  by  our  1979  model.  Based  on  a  constant  release 
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TABLE  3-2 


Change  in  AER  1-0  Model  Calculation  since 
1979  of  Calculated  Effect  on  Total  Ozone 
from  Constant  Release  of  CFC-11  and  CFC-12* 


Model 


Change  in  03  (%) 
Present  Day  Steady  State 


A.  NASA  (1979)  Chem 


-  1.7  -  16.4** 


B.  Same  as  A  with  -  1.6  -  16.0 

k(0H  +  H202)  «  2. 5  (-12) exp 


C.  Same  as  B  with  -  -  11.0 

k(0H  +  HN03)  -  1. 5  (-14)  exp  (£|2) 


D-  NASA  (1980)  Chem  -  -  8.2 

(JPL  Publ.  81-3) 


E.  Same  as  D  with  -  -  4.7 

k(0H  +  HN04)  «  4  (-12) 

k(0H  +  H02)  =  8  (-11) 

F.  NASA/WMO  (1981)  Chem  -  0.65  -  6.1 


*  The  constant  release  rates  of  CFC-11  and  CFC-12  are 
set  at  1977  production  levels. 

**  All  calculations  (except  for  model  F)  use  fast  rate 
for  CINO^  formation. 
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CHANGE  IN  LOCAL  OZONE  ( 10n  molecules  cm*5) 


Figure  3-1 

Comparison  of  1-D  calculated  changes  in  local 
ozone  due  to  CFCs'  perturbation  at  steady  state 
using  different  reaction  rate  constants  from 
NASA  (1979)  and  NASA/WMO  (1981). 
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1960 


YEAR 


1990 


2000 


Figure  3-2 

Comparison  of  1-D  model  calculated  time-dependent  column 
ozone  depletions  due  to  CFC's  release  using  reaction  rate 
constants  from  NASA  (1979)  and  NASA/WMO  (1981). 
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scenario  of  CFCs  at  1979  release  rates  after  1980,  our  model 
with  1981  rates  projects  that  the  column  03  reduction  will 
reach  1.7%  by  year  2000. 

Figures  3-3  and  3-4  show  the  possible  changes  in  03 
caused  by  a  doubling  of  N20  and  CC>2  respectively.  In  con¬ 
trast  to  CFCs  and  N20,  an  increase  in  C02  is  believed  to 
increase  column  03  (Luther  et  al.,  1977;  Boughner,  1978; 
Groves  et  al.,  1978)  as  the  C02-induced  cooling  effect  re¬ 
duces  the  efficiency  of  several  0^  removal  processes.  The 
quantitative  effect  of  C02  on  depends  on  the  absorption 
band  treatment  used  in  the  radiative  calculation  as  well  as 
the  treatment  of  hydrostatic  adjustment.  The  narrow  band 
model  predicts  a  smaller  increase  of  column  03  (2.5%)  than 
that  (3.9%)  predicted  by  the  wide  band  model  due  to  C02 
doubling  (Wang  et  al.,  1982).  It  is  interesting  to  note 
that  the  maximum  C02  impact  on  03  (Figure  3-4)  is  predicted 
to  occur  at  around  40  km,  a  region  where  CFCs  impact  on  03 
is  also  found  to  be  the  greatest  though  opposite  in  sign. 

Ozone  may  also  be  perturbed  by  the  emission  of  NO  into 
the  upper  troposphere  by  subsonic  aircraft.  Injection  of 
NOx  into  this  region  (9-11  km)  may  promote  ozone  production 
through  processes  such  as: 


of  1-D  model  calculated  changes  in  03  due  to  doubling  of  N30 
ion  rate  constants  from  NASA  (1979)  and  NASA/WMO  (1981) .  Z 
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WIDE  BAND 
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Figure  3-4 

Effect  of  absorption  band  treatments  on  the  one-dimensional 
radiative-convective-photochemical  model  calculated  03  change 
due  to  CO2  increase  from  330  to  660  ppmv.  The  background  C1X 
abundance  is  assumed  to  be  1  ppbv.  The  reaction  rate  constants 
are  from  NASA/WMO  (1981)  . 
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Model  calculations  by  the  Lawrence  Livermore  National 
Laboratory  (LLNL)  indicates  that  the  current  operation  of 
subsonic  aircraft  may  contribute  -  about  1%  increase  in  column 
O3.  Much  of  the  calculated  O3  increase,  however,  occurs  in 
the  upper  troposphere,  a  region  where  uncertainties  in  trans¬ 
port  and  photochemistry  are  particularly  large. 

Increasefin  atmospheric  CH^  might  also  perturb  O3  and 
stratospheric  chemistry.  Rasmussen  and  Khalil  (1981) ,  sug¬ 
gest  an  annual  increase  in  CH4  of  about  2ft  per  year  on  the 
basis  of  22  months  of  almost  continuous  measurements  of  CH4 
at  Cape  Meares.  Should  this  trent  continue  in  the  future, 
an  03  increase  on  the  order  of  a  few  percent  is  possible 
(Wang  et  al.,  1980)  since  CH4  oxidation  is  a  source  of  O3 
(Johnston  and  Podolski,  1978)  and  its  reaction  with  Cl  plays 
a  dominant  role  in  the  remrval  of  stratospheric  Cl. 
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3.2  Multiple  Perturbation  Studies 

It  is  clear,  from  the  above  discussions,  that  human 
influences  on  global  ozone  are  diverse  and  often  complex. 

It  seems,  therefore,  unrealistic  to  assess  the  anthropo¬ 
genic  effect  on  ozone  based  on  a  single  perturbation.  This 
is  particularly  true  when  the  interaction  of  perturbations 
is  nonlinear  and  some times  even  opposite  in  sign.  Multiple 
perturbation  studies  have  been  initiated  by  Logan  et  al. 
(1978)  and  Sze  et  al.  (1978).  Their  work  focuses  on  the 
coupled  effect  of  N20  and  CFCs  on  03*  They  note  that  the 
effect  of  individual  perturbations  is  not  simply  additive. 

For  instance,  the  impact  of  CFCs  on  ozone  strongly  depends 
on  the  atmospheric  level  of  N20.  At  an  elevated  N20  level, 
the  calculated  response  of  O3  to  CFC  releases  diminishes 
significantly. 

Table  3-3  illustrates  some  of  the  examples  of  multiple 
scenarios  calculated  by  the  AER  one-dimensional  model  using 
NASA  (1981)  rates.  In  the  joint  scenario  (A)  involving 
simultaneous  doubling  of  N2O  and  constant  emission  rates  of 
CFCs,  the  net  total  ozone  reduction  is  13%,  only  somewhat 
higher  than  the  value  of  10%  predicted  for  doubling  N20 
alone.  It  should  be  noted,  however,  that  the  evidence  of 
an  increasing  trend  in  N20  is  mainly  based  on  Weiss's  (1981) 
observations.  A  long  time  series  is  needed  to  quantify  the 
rate  of  N20  increase  in  the  atmosphere. 
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Model 


Change  in  Column  03  (%) 


A. 

CFCS 

(steady  state) 

+  2  x  N20 

-  13.0 

B. 

CFCS 

(steady  State) 

+  2  X  C02 

-  4.8** 

C. 

CFCS 

(2010)  +  1.5  x 

co2 

-  1.3** 

*  Refer  to  Table  3-1  for  baseline  calculations. 

**  Based  on  narrow  band  model  in  the  C02  radiative 
cooling  calculation. 
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Atmospheric  C02,  on  the  other  hand#  is  observed  to  be  * 
increasing  steadily  since  1958  (Keeling  et  al.,  1976a, b) . 

A  doubling  in  C02  could  occur  around  the  year  2025  to  2075 
(CEQ,  1980) ,  about  100-150  years  before  CFCs  are  projected 
to  attain  their  steady  state  concentration.  The  joint  ef¬ 
fect  of  doubling  C02  and  constant  emission  rates  of  CFCs  on 
total  ozone  is  a  4.8%  reduction,  compared  to  6.1%  calculated 
with  CFCs  emission  alone. 

Scenario  C  represents  a  projection  of  coupled  effect  of 
C02  and  CFCs  on  ozone  on  an  intermediate  time  scale  ('v  year 
2010) .  This  scenario  includes  500  ppm  C02  and  4  ppb  C1X 
[approximately  75%  (i.e.,  3  ppb)  comes  from  CFCs]  in  the 
model.  The  net  reduction  in  total  ozone  is  1.3%,  almost  a 
factor  of  two  smaller  than  the  value  2.3%  calculated  with 
no  C02  increase  (fix  C02  at  current  level  of  330  ppm) . 

3.3  2-D  Perturbations 

The  results  presented  in  Sections  3.1  and  3.2  refer  to 
an  average  global  perturbation  to  calculated  by  1-D  models 
Earlier  calculations  by  Pyle  (1978)  suggest  that  the  maximum 
depletion  of  03  by  CFCs  occurs  in  winter  at  high  latitudes, 
where  the  transmission  of  ultraviolet  radiation  is  negligible 
In  contrast,  the  percentage  of  03  decrease  is  minimal  in 
summer  at  low  latitudes  where  the  penetration  of  ultraviolet 
radiation  is  greatest.  Since  the  exposure  of  human  popula¬ 
tion  and  many  plants  to  ultraviolet  radiation  is  likely  to 
be  seasonal  and  latitudinal  dependent,  a  global  average  03 
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reduction  calculated  by  1-D  models  (including  those  given 
in  Sections  3.1  and  3.2)  may  not  be  very  meaningful,  partic¬ 
ularly  when  2-D  model  calculation  suggests  that  O3  depletion 
may  vary  significantly  with  latitudes  and  seasons. 

Figure  3-5  shows  the  present  day  0^  depletion  by  CFC's 
release  calculated  by  the  AER  2-D  model  with  NASA/WMO  1981 
rates.  A  detailed  description  of  the  model  is  given  in 
Appendix  A.  The  calculated  latitudinal  and  seasonal  depen¬ 
dence  of  depletion  is  qualitatively  similar  to  the  earlier 
calculation  by  Pyle  (1978) .  The  largest  0^  depletion  is  found 
in  the  winter  at  high  latitudes.  Our  calculation,  however, 
differs  quantitatively  from  that  of  Pyle  (1978),  in  one  import¬ 
ant  aspect:  that  the  ratio  of  polar  to  tropical  03  depletion 
is  about  1.8,  a  value  significantly  less  than  the  correspon¬ 
ding  ration  of  3.5  calculated  by  Pyle  (1978).  We  believe 
the  principal  cause  for  the  difference  arises  from  the  much 
faster  rates  for  OH  +  HNO3  and  OH  +  HN04  used  in  our  present 
model.  Because  of  the  higher  stratospheric  concentrations 
of  HNO^  and  HNOj  at  high  latitudes,  these  two  reactions  tend 
to  suppress  more  OH  at  high  latitudes  than  at  the  tropics. 
Reduced  OH  implies  a  less  efficient  chlorine  cycle  (1-2)  as 
discussed  in  Section  2. 

Figure  3-6  shows  the  possible  effect  of  doubling  NjO 
on  column  ozone  calculated  by  our  2-D  model.  The  calculated 
column  03  depletion  shows  large  latitudinal  gradients,  rang¬ 
ing  from  about  5-7%  reduction  in  the  tropics  to  about  25%  in 
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the  polar  regions.  The  calculated  seasonal  gradient  of  0-j 
depletion,  however,  is  relatively  small.  It  should,  however, 
be  pointed  out  that  since  our  model  calculates  excessive  N02 
at  high  latitudes  during  winter  season  (possible  due  to  the 
use  of  24-hour  average  photolysis  rates) ,  it  is  likely  that 
the  ozone  reduction  by  N0X  is  somewhat  overestimated  above 
30°  latitude. 

In  summary,  potential  ozone  perturbations  by  human 
activities  may  not  be  uniform  in  space  and  time.  Some  care 
is  therefore  needed  in  the  interpretation  of  the  deple¬ 
tion  calculated  by  1-D  models,  particularly  in  regards  to 
the  attempt  of  relating  increase  in  UV-B  dose  with  global 
03  depletion. 
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4.  TWO-DIMENSIONAL  MODEL  STUDIES 


4.1  Introduction 

The  AER  zonal-mean  model  is  designed  to  simulate  the 
zonal-mean  behavior  of  atmospheric  trace  gases.  The  domain 
of  the  model  covers  the  globe  from  pole  to  pole  and  up  to  an 
altitude  of  60  km.  The  horizontal  and  vertical  resolutions 

3 

are  approximately  10  Ion  and  3.5  km,  respectively.  The 
dynamic  variables  are  governed  by  the  system  of  equations 
obtained  by  taking  the  Eulerian  zonal-mean  of  the  simplified 
primitive  equations  which  assumes  geotrophic  wind  balance. 
Dynamic  transport  of  trace  gases  is  affected  by  advection 
from  zonal-mean  wind  and  eddy  transport  parameterized  by  a 
symmetric  diffusion  tensor.  The  photochemical  scheme  used 
in  the  2-D  model  is  based  on  the  grouping  techniques  and  cal¬ 
culates  over  35  species  affecting  the  balance.  A  24-hour 
averaging  of  the  photolysis  rate  is  used  to  simulate  diurnal 
average  conditions.  A  more  detailed  description  of  the  model 
is  given  in  Appendix  A.  In  Sections  4.2  and  4.3  the  results 
of  the  calculation  from  this  model  are  presented. 

In  addition  to  this  basic  model  which  was  developed  over 
the  past  3  years,  we  have  implemented  additional  features 
over  the  past  contract  year.  These  include  a  diurnal  package 
for  calculating  the  diurnal  behavior  of  the  short-lived  spe¬ 
cies  and  a  semi-implicit  scheme  for  solving  the  trace-gas 
equation.  These  improvements  together  with  some  anticipated 
changes  will  be  discussed  in  Section  4.4. 
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4.2  Comparison  of  2-D  Models  with  Observations 

Considerable  advances  have  occurred  in  the  area  of 
measurement  of  stratospheric  trace  gases.  The  developments 
have  provided  a  wealth  of  new  data  on  the  concentration  of 
many  important  stratospheric  gases  including  the  long-lived 
species  (N20,  CH4,  CH-jCl,  C2Hg ,  CFCl^,  CF2C12)  of  tropo¬ 
spheric  origin,  the  hydrogen  radicals  H0X  (OH,  H02 ,  H202) 
derived  *rom  the  photodecomposition  of  stratospheric  H20, 
the  nitrogen  radicals  N0X  (NO,  N02,  HNO^)  from  the  decompo¬ 
sition  of  N20  and  the  chlorine  radicals  (Cl,  CIO,  HC1,  ClNO-j) 
from  the  photodecomposition  of  a  variety  of  halocarbons. 

In  the  area  of  stratospheric  modeling,  considerable 
progress  has  also  been  made  in  the  development  of  a  two- 
dimensional  (2-D)  model,  although  the  assessment  of  potential 
environmental  effects  from  human  activities  still  relies  mainly 
on  one-dimensional  (1-D)  models  (NASA,  1979?  NAS,  1979). 

The  major  strength  of  1-D  models  lies  in  their  ability  to 
include  a  detailed  description  of  chemical  processes  without 
excessive  demands  on  computer  time.  The  problematic  aspects 
of  1-D  models  concerns  the  interpretation  of  their  simulated 
profiles  which  are  considered  to  represent  some  sort  of  global 
(longitudinal  and  latitudinal)  averages.  Strictly  speaking, 
1-D  model  results  should  be  compared  with  the  observed  global 
averages  defined  by  a  large  set  of  spatial  and  temporal  data. 
Acquisition  of  such  a  data  set  for  any  trace  gas  is  clearly 
an  enormous  task  -  a  task  unlikely  to  be  completed  in  the 
foreseeable  future. 


41 


The  main  advantage  of  the  2-D  model  lies  in  its  capa¬ 
bility  of  simulating  atmospheric  latitudinal  gradients  and 
seasonal  behavior  by  accounting  for  the  important  time-depen¬ 
dent  meridional  transport  in  the  atmosphere,  together  with 
the  seasonal  and  latitudinal  variations  of  the  solar  zenith 
angle.  The  recent  improvement  in  2-D  models  involves  the 
inclusion  of  detailed  chemical  processes.  For  instance,  the 
chemical  content  of  DuPont's  (Miller  et  al.,  1981)  and  AER's 
(Ko  et  al.,  1980)  2-D  model  is  now  at  par  with  that  of  the 
most  advanced  1-D  models. 

On  the  other  hand,  one  should  recognize  the  fact  that 
current  2-D  models  are  not  developed  from  first  principles. 
Rather,  the  transport  in  the  models  is  partly  parameterized 
by  eddy  diffusion  coefficients  derived  from  observations. 

It  is  therefore  essential  to  substantiate  the  models  by 
comparison  with  observations.  Specifically,  observations 
complement  model  results  in  at  least  four  general  areas: 
i)  boundary  conditions,  ii)  atmospheric  transport,  iii) 
chemical  system  and  iv)  combined  chemistry  and  transport. 

First,  observations  often  are  used  to  define  the  nec¬ 
essary  boundary  conditions  posed  in  the  2-D  models.  Second, 
because  of  the  simplicity  of  chemistry  for  certain  long-lived 
species  such  as  NjO,  CH^  and  CFCs,  data  on  these  gases  may  be 
used  to  check  on  the  appropriateness  of  transport  parameter¬ 
ization.  Third,  data  on  radical  species  may  be  used  to  test 
our  current  understanding  of  various  chemical  cycles  since 
the  partition  among  the  radical  species  is  generally  insen- 
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sitive  to  transport  although  transport  may  affect  the  abso¬ 
lute  abundance  of  the  radical  calculated  in  the  model. 
Finally,  the  observed  spatial  and  temporal  ozone  distribu¬ 
tion,  which  is  believed  to  be  controlled  by  both  chemistry 
and  transport,  may  be  used  to  test  our  overall  understanding 
of  the  stratospheric  chemical  and  dynamical  system.  The  last 
step  obviously  depends,  to  a  large  extent,  on  how  well  the 
model  results  agree  with  the  observed  long-lived  and  radical 
species. 

This  section  is  organized  in  the  manner  of  discussing 
each  of  the  above  areas  of  interst.  Theoretical  results 
presented  in  this  section  are  based  on  the  AER  2-D  model 
with  NASA/WMO  1981  rates. 

Our  numerical  experiment  was  carried  out  in  two  stages. 
We  first  performed  the  calculation  solving  for  the  annual 
periodic  steady  state  distributions  of  ^O,  CH4,  I^,  CH^Cl, 
CC14,  CH3CC14,  CO,  03,  Clx  and  N0X  using  fixed  mixing  ratio 
conditions  for  the  lower  boundary  as  summarized  in  Table  A-4 
in  Appendix  A.  The  resulting  concentrations  can  be  inter¬ 
preted  as  those  of  the  "clean"  atmosphere  prior  to  introduc¬ 
tion  of  CFC-11  and  CFC-12.  The  model  atmosphere  is  then 
integrated  in  time  for  30  model  years  where  CFC-11  and  CFC- 
12  are  introduced  into  the  model  atmosphere  with  time-depen- 
dent  flux  boundary  conditions  using  the  em^-jion  rates  for 
the  years  1950  to  1980  (CMA,  1980).  The  calculated  values 
correspond  to  the  gas  concentration  of  the  present  day  atmo¬ 
sphere. 
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4 . 2a  Long-Lived  Species 

Among  the  long-lived  species,  the  upward  diffusing 
species  -  N20,  CH4,  CC14,  CH3CI,  CH-jCCl-j,  CFC-11  and  CFC-12  - 
are  characterized  by  the  fact  that  there  is  no  in  situ  pro¬ 
duction  in  the  atmosphere  for  these  molecules.  The  gases 
enter  the  stratosphere  by  upward  diffusion  from  the  tropo¬ 
sphere  and  are  removed  by  photolysis,  0(1D)  and  OH  reactions 
in  the  stratosphere.  An  observed  feature  common  among  the 
upward  diffusing  species  is  that  at  each  altitude  above  the 

tropopause,  higher  concentrations  are  always  found  in  the 
i-r°pical  stratosphere.  This  is  to  be  expected  from  the  known 
meteorological  patterns  through  which  upward  transport  occurs, 
preferentially  in  the  tropical  latitudes,  followed  by  poleward 
transport  within  the  stratosphere.  This  important  feature  is 
present  in  our  model  calculation  for  the  aforementioned  gases 
and  is  illustrated  by  a  latitude-altitude  plot  of  the  mixing 
ratio  for  N20  in  Figure  4-1.  The  figure  shown  corresponds  to 
winter  conditions  in  the  northern  hemisphere  (December) .  It 
is  interesting  to  note  that  there  is  a  poleward  rise  in  the 
stratospheric  concentration  in  the  winter  hemisphere  at  high 
lati-t*ides  reflecting  the  lack  of  chemical  removal.  This  model 
feature,  of  course,  remains  to  be  verified  by  future  observa¬ 
tions. 

Among  the  upward  diffusing  species,  sufficient  data  from 
different  latitudes  are  available  for  N20,  CFC-11  and  CFC-12. 
The  bulk  of  these  data  was  collected  by  four  research  groups: 
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Figure  4-1 

titude-altitude  cross-sections  for  the  calculated 
xing  ratio  of  NjO  for  the  month  of  December. 
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NOAA  Goldan  et  al.  (1980,  1981) .  Balloon-borne 
evacuated  grab  samples. 

NCAR  Heidt  et  al.  (1975) .  Balloon-borne  cryosampler. 

KFA  Schmidt  et  al.  (1981).  Balloon-borne  cryosampler. 

Ames  Tyson  et  al.  (1978);  Vedder  et  al.  (1978);  Vedder 
et  al.  (1981) .  Balloon-  and  aircraft-borne  cryo¬ 
sampler. 

n2° 

Figures  4-2a  and  b  show  the  calculated  NjO  profiles 
together  with  observations  for  the  equatorial  and  mid-lati¬ 
tude  regions  respectively.  Note  that  the  data  are  collected 
over  a  number  of  flights  and  this  may  account  for  the  scatter 
in  the  data  at  high  altitudes.  The  data  show  that,  at  any 
latitude,  the  mixing  ratios  of  N20  decrease  with  altitude 
above  the  tropopause;  a  feature  reflecting  that  these  gases 
are  removed  in  the  stratosphere.  Another  interesting  feature 
concerns  the  latitudinal  variations  of  these  species  at  con¬ 
stant  pressure  (altitude)  levels.  Specifically,  the  data 
show  that  the  mixing  ratios  at  a  given  altitude  (above  tro¬ 
popause)  decrease  towards  the  high  latitude  regions. 

The  agreement  between  the  2-D  model  calculations  and 
observations  for  N^O  is  good.  The  model  successfully  simu¬ 
lates  the  observed  vertical  gradients  as  well  as  the  lati¬ 
tudinal  gradients  of  N2O.  As  discussed  in  Appendix  A,  the 
dynamic  transport  in  the  AER  model  uses  the  advection  wind- 
fields  calculated  by  Harwood  and  Pyle  (1977)  and  the  eddy 
diffusion  coefficients  from  Luther  (1973) .  No  attempt  was 
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Ml  PLATITUDES 


Comparison  of  calculated  N2O  profiles  with  observations.  The 
data  are  from  the  measurement  programs  of  NOAA,  NCAR  and  NASA 
Ames  as  compiled  for  the  1981  NASA  Workshop  publication. 


made  to  adjust  any  of  the  parameters.  Thus,  the  favorable 
agreement  between  observation  and  calculation  should  perhaps 
be  viewed  as  a  partial  validation  of  the  dynamic  parameters 
used. 


CFC-11  and  CFC-12 

Since  1975,  a  large  number  of  measurements  of  CFC-11 
and  CFC-12  have  been  reported.  Most  measurements  involve 
laboratory  analysis  by  gas  chromatography  of  air  samples 
collected  by  balloon-borne  samples.  The  research  groups 
which  have  produced  large  numbers  of  data  points  are  the 
four  groups  mentioned  earlier  in  the  section,  viz.  NOAA, 
NCAR,  KFA  and  NASA  Ames.  The  NASA  Workshop  (June,  1981) 
has  recommended  that  some  of  the  NCAR  data  be  excluded 
because  of  problems  associated  with  the  sampling  procedure. 
In  the  mid-latitude  region,  only  those  from  NOAA  and  KFA 
will  be  considered. 

The  data  were  collected  over  a  period  of  time  between 
1975  and  1979.  The  concentrations  of  CFC-11  and  CFC-12 
change  with  time  since  the  gases  are  introduced  with  time- 
dependent  input  rates.  Figure  4-3  shows  the  normalized 
concentrations  for  CFC-11  and  CFC-12  as  functions  of  time 
for  different  altitudes.  Note  that  the  concentrations  at 
25  km  increase  by  about  a  factor  of  1.6  from  1975  to  1979. 
Thus,  the  data  should  be  compared  with  the  calculations  for 
the  appropriate  time. 
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1  calculation  concentrations  for  CFC-11  and  CFC-12  as  a  function 
ime.  The  concentrations  for  2  levels  are  plotted;  each  of  them 
normalized  with  the  1970  concentrations  equal  to  1.  The  values 
taken  from  the  2-D  model  calculation  at  40°  N. 
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Data  for  CFC-11  and  CFC-12  are  plotted  in  Figures  4-4 
through  4-7  together  with  the  calcualted  profiles.  The 
altitude  profiles  are  from  the  time-dependent  calculation 
corresponding  to  the  years  1975  and  1979  as  indicated. 
Agreement  is  reasonably  good  for  CFC-11  at  the  equator  but 
poor  at  mid-latitude.  The  calculated  mid-latitude  CFC-11 

mixing  ratios  above  25  km  lie  significantly  above  NOAA's 
data.  Although  the  overall  agreement  is  somewhat  better 
for  CFC-12,  the  calculated  mid-latitude  CFC-12  mixing  ratios 
above  25  km  are  higher  than  observations.  In  sum,  the  model 
seems  to  have  underestimated  the  latitudinal  gradients  for 
CFC-12  and  more  so  for  CFC-11. 

There  are  several  possible  explanations  for  the  poor 
agreement  between  model  and  observation.  First,  the  photo- 
lytic  destruction  rates  for  CFC-11  and  CFC-12  may  be  under¬ 
estimated  in  the  model  because  of  the  errors  in  the  treat¬ 
ment  of  the  Schumann  Runge  bands.  Second,  the  transport 
rates  may  be  in  error.  Specifically,  the  diffusion  coeffi¬ 
cient  K  which  controls  the  latitudinal  gradients  of  the 
mixing  ratios  may  be  too  fast.  A  slower  K  may  increase 
the  magnitude  of  the  latitudinal  gradient. 

Since  CFC-12  and  N20  are  chemically  more  inert  than 
CFC-11  in  the  stratosphere,  we  therefore  expect  that  the 
calculated  stratospheric  distribution  of  CFC-11  is  more 
sensitive  to  the  adjustment  of  diffusion  coefficients.  In 
other  words,  it  is  possible  to  tune  diffusion  coefficients 
for  better  agreement  between  calculated  and  observed  CFC-11, 
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Figure  4-4 


Comparison  of  the  calculated  CFCI3  profiles  with 
observations.  The  data  are  from  the  measurement 
programs  of  NOAA,  NCAR  and  NASA  Ames  as  compiled 
for  the  1981  NASA  Workshop  publication. 

-  Equator  -  corresponding  to  1979  concentration 

-  -  -  Equator  -  corresponding  to  1975  concentration 
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Figure  4-5 

Comparison  of  the  calculated  CFCI3  profiles  with 
observation.  The  data  are  from  the  measurement 
programs  of  KFA  and  NOAA  as  compiled  for  the  1981 

•  NASA  Workshop  publication.  The  profiles  ccrres- 

pond  to  1975  and  1979  concentrations  in  Figures  a 
and  b  respectively. 

-  30°  N 

- 40°  N 

•  50°  N 
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Figure  4-6 

Comparison  of  the  calculated  CF2CI0  profiles  with 
observations.  The  data  are  from  the  measurement 
programs  of  NOAA,  NCAR  and  NASA  Ames  as  compiled 
for  the  1981  NASA  Workshop  publication. 

-  Equator  -  corresponding  to  1979  concentration 

-  -  -  Equator  -  corresponding  to  1975  concentration 


Figure  4-7 

Comparison  of  the  calculation  of  CF2Cl2  profiles 
with  observations.  The  profiles  correspond  to 
1975  and  1979  concentrations  in  Figures  a  and  b 
respectively. 
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with  only  relatively  small  changes  in  the  calculated  profiltes 
of  NpO  and  CFC-12. 

CH4  and  CH3C1 

Both  CH^  and  CH^Cl  differ  from  the  species  discussed 
previously  in  that  they  are  removed  in  the  troposphere  by 
reaction  with  OH.  Reactions  with  0(^D)  and  Cl  also  play  a 
role  in  the  stratospheric  removal  of  CH4 . 

The  agreement  between  the  calculated  CH4  distribution 
(Fig.  4-8)  and  observation  is  fairly  good,  although  the  cal¬ 
culated  mixing  ratio  lies  somewhat  below  the  observations 
in  the  equatorial  region.  It  should  be  noted  that  the  data 
are  based  on  a  single  observation  profile.  A  smaller  strato¬ 
spheric  OH  abundance  or  a  more  efficient  upwelling,  via  a 
faster  K  ,  for  example,  could  bring  the  model  into  closer 

«  A 

agreement  with  observations. 

The  calculated  mid-latitude  profiles  for  CH^Cl  (Fig.  4-9) 
lie  above  the  observations.  It  must  be  pointed  out  that  there 
are  less  data  available  for  CH^Cl  than  for  CFC-11  and  CFC-12. 
Nevertheless,  the  behavior  of  the  model  calculated  CH^Cl 
appears  to  be  consistent  with  those  of  CFC-11  and  CFC-12 
pointing  to  the  possibility  that  a  slower  K  may  increase 
the  latitudinal  gradient  and  improve  the  agreement  with 
observations. 
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Figure  4-8 

Comparison  of  the  calculated  CH4  profiles  with 
observations  for  different  latitudes.  The  data 
are  from  Bush  et  al.  (1979),  Heidt  and  Krasnec 
(1978),  Ehhalt  (1978)  and  Fabian  et  al.  (1979) 
as  compiled  in  NASA  1049.  The  curves  (-  -  -, 

- )  correspond  to  (Equator,  10°  N) ,  (20°  N, 

40°  N)  and  (40°  N,  60°  N)  respecitvely  in  the 
three  panels. 
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Figure  4-9 

Comparison  of  the  calculated  CH3CI  profiles  with 
observations  for  different  latitudes.  The  data 
from  Penkett  et  al.  (1980)  and  Schmidt  et  al. 
(1960)  are  compiled  for  the  1981  NASA  Workshop 
The  curves  (-.-.-)  30°  N,  (-  -  -)  40°  N  and 
( - )  50°  N. 


H2  and  CO 


The  other  long-lived  species  for  which  there  are  avail¬ 
able  data  are  H2  and  CO.  The  interpretation  of  the  compari¬ 
son  between  calculations  and  observations  is  less  straight 
forward  because  these  species  are  also  produced  in  the  atmo¬ 
sphere  via  photochemical  reactions. 

For  H2*  the  agreement  between  observation  and  calcula¬ 
tion  is  good  below  30  km  (Fig.  4-10) .  Above  30  km#  the  cal¬ 
culated  profiles  lie  below  the  observations.  It  should  be 
noted  that  the  observations  show  less  of  a  latitudinal  gra¬ 
dient  than  the  upward  diffusive  species, possibly  indicating 
that  the  effect  of  in  situ  chemical  production  may  be  com¬ 
parable  to  dynamic  transport  in  determining  local  concentra¬ 
tions  . 

The  amount  of  CO  data  is  rather  limited.  Figure  4-11 
shows  the  calculation  together  with  data.  The  agreement 
is  good  except  for  the  data  above  30  km  which  are  based  on 
column  measurements  of  Zander  et  al.  (1981) . 

4 . 2b  Radical  Species 

The  radical  species  are  usually  characterized  by  chem¬ 
ical  lifetimes  much  shorter  than  the  transport  time  scale. 

The  partition  among  various  radical  species  is  therefore 
mainly  controlled  by  photochemistry,  although  their  abso¬ 
lute  concentrations  still  may  be  affected  by  the  stratospheric 
distribution  of  the  long-lived  (precursor)  species  discussed 
in  the  earlier  section.  Nevertheless,  a  comparison  between 
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Figure  4-10 

Comparison  of  the  calculated  H2  profiles 
with  observations.  The  data  are  compiled 
for  the  1981  NASA  Workshop  publication. 
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the  observed  and  calculated  radicals  at  different  altitudes, 
latitudes  and  seasons  should  provide  useful  insights  into 
the  stratospheric  chemical  system,  particularly  when  the  dis¬ 
tribution  of  the  long-lived  species  is  simulated  well  in  the 
model. 

The  vertical  profiles  for  different  latitudes  and  sea¬ 
sons  presented  here  correspond  to  the  daytime  condition. 
Quantitative  comparison  of  the  calculated  radical  concen¬ 
trations  with  observation  remains  difficult  for  the  follow¬ 
ing  reasons.  First,  the  concentrations  are  often  very  sensi¬ 
tive  to  the  detailed  diurnal  variation  of  solar  zenith  angle. 
A  full  diurnal  treatment  is  not  included  in  our  2-D  calcula¬ 
tion  because  of  the  enormous  computing  requirement.  Second, 
observations  may  be  subject  to  local  variabilities  influenced 
by  factors  not  included  by  current  models.  The  calculated 
latitude  and  altitude  cross-section  of  the  radical  species 
are  given  in  Appendix  B.  Note  that  those  calculations  were 
performed  using  24-hour  average  conditions.  The  result  pre¬ 
sented  in  Appendix  B  is  useful  in  giving  a  perspective  on 
the  latitudinal  variation  of  the  species.  For  comparison 
with  data,  it  is  more  convenient  to  use  altitude  profiles. 

A  common  feature  among  the  calculated  profiles  is  that  the 
summer  profiles  for  different  latitudes  are  much  closer  to¬ 
gether  than  the  corresponding  set  of  winter  profiles  in  the 
northern  hemisphere.  This  feature  appears  to  reflect  the 

smaller  differences  in  the  solar  zenith  angle  among  the  dif¬ 
ferent  latitudes  in  the  summer  hemisphere. 
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03(P),  OH,  H02 

Figures  4-12,  4-13  and  4-14  show  the  calculated  profiles 
for  03(P),  OH  and  H02  respectively  together  with  the  measure- 

3 

ments.  The  agreement  for  0  (P)  is  good.  The  OH  and  H02 
profiles  are  calculated  using  a  uniform  concentration  of 
5  ppmv  for  stratospheric  water  vapor.  The  agreement  for  OH 
is  good.  In  the  case  of  HO 2,  the  calculated  profiles  are 
somewhat  lower  than  the  observations.  Note  that  the  data 
for  OH  and  H02  below  30  ]cm  are  still  non-existent. 

NO,  N02,  HN03 

The  calculated  behavior  of  the  nitrogen  radicals  is 
affected  by  NOx  distribution.  Figure  4-15  shows  the  calcu¬ 
lated  NOx  profiles  for  different  latitudes  corresponding  to 
summer  condition. 

Our  calculations  show  that  in  the  upper  stratosphere, 
the  bulk  of  the  NOx  is  in  the  form  of  NO.  Note  that  the 
latitudinal  variation  of  NO  (Fig.  4-16)  closely  mimics  that 
of  NOx  in  that  region.  In  the  lower  stratosphere,  most  of 
the  NOx  is  in  the  form  of  N02  and  HNO3.  The  behavior  of  the 
N02  profile  for  50°  N  in  winter  (Fig.  4-17)  seems  contrary 
to  the  recent  seasonal  data  reported  by  Noxon  (1979)  and 
Coffey  et  al.  (1981)  which  showed  that  there  should  be  less 

N02  in  the  winter  hemisphere  at  high  latitudes.  Our  prelim¬ 
inary  investigation  shows  that  this  may  have  resulted  from 
the  diurnal  treatment  being  used.  A  more  realistic  treat- 
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Figure  4-15 

Calculated  altitude  profiles  for  N0X  and  N2O  for  different 
latitudes  corresponding  to  summer  in  the  northen  hemisphere. 
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Calculated  day  time  concentrations  of  NO2  in  northern  hemisphere  for  different 
seasons.  The  shaded  areas  correspond  to  data  reported  by  Ackerman  and  Muller 
(1973),  Murcray  et  al.  (1974),  Ackerman  et  al.  (1975),  Fontanella  et  al.  (1975), 
Kerr  and  McElroy  (1976),  Evans  et  al.  (1978),  Goldman  et  al.  (1978)  and  Blather- 
wick  et  al.  (1980) . 


ment  of  the  diurnal  behavior  may  be  required  to  simulate  the 
high  latitude  NC>2  during  winter  conditions. 

The  calculated  profiles  of  HNO^  agree  well  with  obser¬ 
vation  in  the  lower  stratosphere  but  lie  above  observations 
in  the  upper  stratosphere,  as  shown  in  Figure  4-18.  Obser¬ 
vations  of  HNO-j  column  as  a  function  of  latitude  are  also 
available  and  show  a  minimum  at  the  equator.  Figure  4-19 
shows  the  calculated  HNO^  column  abundance  together  with  the 
observations.  Both  calculations  and  observations  indicate 
that  the  minimum  HNO^  column  occurs  around  the  equator.  The 
behavior  of  column  HNO^  in  Figure  4-19  is  consistent  with 
the  latitudinal  variation  of  NO  in  Figure  4-15. 

On  balance,  the  calculated  concentrations  agree  reason¬ 
ably  well  with  the  observations  except  for  NO  and  HN03  in  the 
upper  stratosphere,  where  the  model  tends  to  overestimate  the 
species  abundance. 


CIO,  HC1,  C1N03,  Cl 

Figure  4-20  shows  the  ClX  profiles  for  different  lati¬ 
tudes  corresponding  to  summer  conditions.  In  contrast  to 
NOx,  ClX  is  well  mixed  in  the  upper  stratosphere  (no  vertical 
gradient)  because  of  the  absence  of  a  stratospheric  sink. 

Figure  4-21a  and  b  show  the  altitude  profiles  for  HCl 
which  is  the  major  stratospheric  chlorine  reservoir.  The 
calculated  HCl  appears  to  lie  within  a  factor  of  2  of  the 
observations,  although  in  the  upper  stratosphere  (^40  km) , 
the  model  results  like  in  the  lower  range  of  observations. 
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VERTICAL  COLUMN  DENSITY  AT  12km  OCT  molecules  cm 


Figure  4-19 


Calculated  column  abundance  of  HN03  corresponding  to 
spring  time  condition.  Shaded  area  corresponds  to 
the  data  reported  by  Murcray  et  al.  (1975) . 
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A  higher  stratospheric  chlorine  content  would  improve  the 
agreement. 

Figure  4-22  shows  the  calculated  profiles  of  CIO  to¬ 
gether  with  the  observations.  With  the  recent  revision  in 
the  reaction  rates  as  discussed  in  Section  2,  the  agreement 
in  the  lower  stratosphere  is  much  better  than  before  (prior 
to  1980).  However,  the  observed  mixing  ratios  appear  to  peak 
at  higher  altitudes  than  our  model  does.  Furthermore,  the 
model  tends  to  underestimate  the  CIO  abundance  by  about  a 
factor  of  2  above  35  km.  The  situation  is  somewhat  disturb¬ 
ing,  since  the  discrepancy  occurs  in  the  region  where  the 
CFCs- induced  ozone  depletion  is  predicted  to  be  the  greatest. 
This  apparent  discrepancy  may  be  caused  by  a  combination  of 
factors,  such  as  underestimation  of  upper  stratospheric  OH 
or  total  chlorine  content  by  the  model.  If  the  upper  strato¬ 
spheric  OH  were  underestimated  by  the  model,  then  it  is  pos¬ 
sible  that  current  models  might  have  understated  the  impact 
of  CFCs  on  ozone.  This  interpretation  may  be  invalid  if  the 
discrepancy  is  actually  due  to  the  underestimate  of  the  strato 
spheric  chlorine  content  by  the  models.  It  is  therefore 
difficult  to  relate  the  present  CIO  data  with  CFCs  impact 
without  knowing  the  cause  for  the  discrepancy  between  obser¬ 
vation  and  model.  A  simultaneous  measurement  of  OH,  CH^, 

CIO  and  HC1  around  38  km  will  help  to  resolve  this  i  ,ue. 

The  calculated  profiles  for  CINO^  are  given  in  Figure 
4-23.  The  concentrations  calculated  with  the  slow  CINO^ 
formation  rate  as  recommended  by  NASA/WMO  (1981)  are  signif- 
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Figure  4-22 

Calculated  day  time  concentrations  of  CIO  in  the  northern  hemisphere  for  different 
seasons.  The  curve  labeled  Anderson  1980  compared  to  the  data  of  Anderson  excluding 
the  July  14  measurement. 
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Figure  4-23 


Calculated  day  time  concentrations  of  CINO3 
northern  hemisphere  for  different  seasons, 
data  are  from  Murcray  et  al.  (1979), 


in  the 
The 


icantly  lower  than  the  observations  of  Murcray  et  al.  (1979) . 
Measurement  of  this  species  may  provide  some  insights  as  to 
which  formation  rate  of  CINO^  under  stratospheric  conditions 
is  more  appropriate. 

The  calculated  profiles  for  Cl  are  given  in  Figure  4-24. 
Again,  as  in  the  case  of  CIO,  the  model  seems  to  have  under¬ 
estimated  the  Cl  abundance  at  around  the  38-40  km  region. 

The  possible  causes  for  the  difference  between  calculated 
and  observed  Cl  at  40  km  may  be  similar  to  that  for  CIO  as 
discussed  earlier. 


4 . 2c  Ozone 

Figures  4-25a  and  b  show  the  photochemical  lifetime  of 
ozone  calculated  via  the  expression 


2  x J2102] 


(4.1) 


as  a  function  of  latitude  and  altitude  for  two  seasons. 

Since  the  stratospheric  dynamic  time  scale  is  typically  a 
few  weeks,  the  figure  also  delineates  regions  whre  photo¬ 
chemistry  or  dynamic  transport  dominate.  Typically,  in  the 
tropics  above  30  km,  03  concentration  is  determined  by  local 

chemical  sources  and  sinks,  whereas  in  the  lower  stratosphere 
(z  <  30  km)  03  is  mainly  controlled  by  transport.  Since 
much  of  the  total  ozone  resides  below  30  km,  the  behavior 
of  column  03  must  be  strongly  influenced  by  atmospheric 


transport. 


Calculated  photochemical  lifetime,  of  ozone  as  a  function  of  latitude  and  altitude 
for  different  seasons.  The  lifetime  is  given  by  the  ratio  [0,3 /2J~ [0,] . 


Figure  4-26  shows  the  observed  behavior  of  column  0 ^ 
as  a  function  of  latitude  and  season  as  compiled  by  Oeutsch 
(1971)  from  ozone  measurements  obtained  from  the  Dobson  net¬ 
work.  The  essential  features  are  the  relatively  low  values 
(^260  Dobson  units)  around  the  equatorial  region  and  the 
occurrence  of  the  spring  time  maxima  (April  and  November) 
at  high  latitudes  in  each  hemisphere. 

Figure  4-27  shows  the  calculated  Oj  distribution  of 
the  present  day  atmosphere  from  our  2-D  model.  Though  the 
calculated  features  are  qualitatively  correct,  the  calcu¬ 
lated  latitudinal  gradient  of  column  ozone  is  somewhat  smal¬ 
ler  than  the  observed  values.  In  addition,  the  northern 
hemispheric  spring  time  maximum  occurs  at  a  later  time  (May 
instead  of  late  March) . 

It  is  well  known  that  transport  is  mainly  responsible 
for  the  latitudinal  and  seasonal  behavior  of  column  03,  We 
therefore  performed  a  numerical  experiment  to  see  if  adjust¬ 
ment  of  the  eddy  diffustion  coefficient  would  lead  to  better 
agreement.  In  the  formulation  of  Reed  and  German  (1965) , 

the  value  o<=K  /K  gives  the  slope  of  the  mixing  path.  A 
yx  yy 

higher  value  for  a  would  imply  a  larger  counter  gradient 
transport  from  the  equatorial  to  the  high  latitude  region. 

In  the  experiment,  we  double  the  value  of  Kyz  in  the  model 
calculation.  In  order  to  preserve  numerical  stability,  the 
value  of  K  has  to  be  multiplied  by  4  at  the  same  time. 

The  resulting  0^  column  distribution  is  shown  in  Figure  4-28 
which  shows  a  larger  contrast  between  the  equatorial  minima 
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Figure  4-26 

Observed  ozone  column  as  a  function  of  latitude  and 
season  as  compiled  by  Deutsch  (1971) . 
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Figure  4-27 

Calculated  column  ozone  as  a  function  of  latitude  and 
season.  The  reaction  rate  constants  are  from  NASA/WMO 


and  the  spring  time  maxima.  A  similar  result  has  been  ob-  . 
tained  by  Harwood  and  Pyle  (1980) . 

We  next  performed  several  numerical  experiments  to 
examine  the  effect  of  chemistry  on  the  calculated  O^.  Fig¬ 
ure  4-29  shows  the  result  from  a  model  run  that  incorporates 
only  oxygen  and  hydrogen  chemistry.  As  expected,  the  magni¬ 
tude  of  the  calculated  0^  column  is  higher  than  the  observa¬ 
tion.  However,  the  ratio  of  the  spring  time  maximum  to  the 
equatorial  value  of  column  03  is  much  closer  to  the  observed 
value.  Note  that  the  spring  time  maximum  also  occurs  at 
approximately  the  right  time. 

Figure  4-30  shows  the  calculated  column  when  nitrogen 
chemistry  is  added.  From  the  comparison  of  Figure  4-29  with 
Figure  4-30,  it  is  evident  that  N0X  has  relatively  little 
impact  on  equatorial  column  03  but  very  large  impact  on  high 
latitude  column  03.  This  observation  leads  us  to  explore 
the  possibility  that  the  present  model  with  a  24-hour  aver¬ 
age  scheme  might  have  overestimated  the  03  removal  at  high 
latitudes  due  to  the  reaction  of  N02  +0.  We  are  currently 
incorporating  detailed  diurnal  treatment  in  our  2-D  model  as 
a  possible  means  to  reconcile  model  results  with  observations. 

Figure  4-31  shows  the  calculated  mixing  ratio  of  03  as 
a  function  of  latitude  and  altitude  for  different  seasons. 

The  raaximun  occurs  at  an  altitude  of  35  km,  and  at  a  lati¬ 
tude  corresponding  roughly  to  the  position  of  the  overhead 


sun  for  each  season. 
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Figure  4-28 

Calculated  column  ozone  from  the  numerical  experiment 
testing  the  sensitivity  of  column  ozone  to  changes  in 
the  eddy  diffusion  coefficient.  In  the  experiment,  K 
and  Kzz  increased  by  a  factor  of  2  and  4  respectively 


OZONE  COLUMN  (Dobson  Units) 

Figure  4-30 

Dzone  using  oxygen-hydrogen-nitrogen 
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Figure  4-31 


Calculated  mixing  ratio  of  0^  for  different  seasons. 
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Figure  4-32  shows  the  observed  cross-sections  of  0^ 
for  different  seasons  as  compiled  by  Deutsch  (1978) .  The 
overall  feature  of  our  calculated  distribution  is  similar 
to  the  observation.  The  calculated  maximum  mixing  ratio  of 
12  ppmv  appears  to  be  somewhat  higher  than  the  observations 
which  shows  a  broad  maximum  of  >8  ppmv. 

4 . 2d  Summary  of  Comparisons 

We  have  compared  the  AER  2-D  model  results  with  obser¬ 
vations  of  important  stratospheric  gases.  In  our  previous 
studies  (Sze  et  al.,  1978,  1980),  we  have  pointed  to  several 
discrepancies  between  the  observed  and  the  calculated  concen¬ 
trations  of  CIO,  C2Hg  and  ratios  of  HN02:N02  and  HF:HC1  in 
the  lower  stratosphere  (20-30  km) .  We  suggested  that  these 
discrepancies  might  be  attributed  to  a  common  cause  -  that 
models  have  overestimated  the  OH  concentration  in  the  20-30 
km  region  where  data  for  OH  are  lacking.  Recent  revisions  in 
the  rate  constants  for  several  key  H0X  reactions  have  signif¬ 
icantly  lowered  the  calculated  stratospheric  OH  abundance  and 
thereby  greatly  improved  the  agreement  between  model  and  ob¬ 
servations  except  for  C2Hg  which  still  lies  significantly 
above  observations. 

Several  important  discrepancies,  however,  still  remain. 
The  discrepancy  concerning  the  calculated  and  observed  CIO 
at  around  40  km  has  actually  been  amplified  when  the  new 
rate  data  were  adopted.  The  discrepancy  in  CIO  is  partic¬ 
ularly  disturbing  since  it  occurs  at  a  region  (^40  km) 
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where  catalysis  of  0^  removal  by  CIO  is  predicted  to  be  the 
greatest.  The  precise  cause  for  this  discrepancy  is  currently 
unknown . 

Another  important  discrepancy  concerns  the  excessive 
amount  of  stratospheric  (25-*30  km)  CFCs  (particularly  CFC- 11) 
calculated  by  the  models.  Again,  the  cause  for  this  discrep¬ 
ancy  is  unknown,  although  it  has  been  suggested  that  photol¬ 
ysis  rates  of  CFCs  in  the  lower  stratosphere  might  have  been 
underestimated  by  models.  Resolution  of  this  discrepancy 
could  have  important  implications  for  CFCs'  lifetime  estim¬ 
ates. 

Agreement  between  the  calculated  and  observed  vertical 
profiles  of  N0X  species  at  mid-latitudes  is  reasonably  good, 
except  for  HNC>3  above  30  km,  a  region  where  the  model  calcu¬ 
lates  significantly  more  HNO3  than  observation.  On  the  other 
hand,  our  2-0  model  did  relatively  poorly  in  the  simulation 
of  the  latitudinal  and  seasonal  behaviors  of  NO2  and  NO,  par¬ 
ticularly  in  winter  time  at  high  latitudes.  This  discrepancy, 
however,  remains  to  be  investigated  with  a  2-D  model  with  full 
diurnal  treatment  in  stratospheric  chemistry  before  one  can 
isolate  the  causes  of  such  a  discrepancy. 

Finally,  the  calculated  column  ozone  distributions  are 
qualitatively  consistent  with  observations,  although  the  model 
seems  to  have  underestimated  the  latitudinal  and  seasonal  gra¬ 
dients  of  column  ozone.  In  view  of  the  several  discrepancies 
between  the  calculated  and  observed  trace  gases  discussed 
above,  it  is  perhaps  unrealistic  to  expect  excellent  agreement 
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between  the  calculated  and  observed  column  ozone.  However, 
if  the  above  discrepancies  in  the  trace  gase&  were  resolved, 
then  any  remaining  discrepancy  in  column  ozone  should  be  con¬ 
sidered  very  serious. 

4.3  Model  Predictions  of  Atmospheric  Lifetimes  for  N~0, 

CFC-11  and  CFC-12  * 

As  evident  from  observation  as  well  as  from  our  model 
results  (Cf.  Fig.  4-2,  4-4,  4-5,  4-6,  4-7),  the  spatial  dis¬ 
tribution  of  ^0,  CFC-11  and  CFC-12  varies  significantly  with 
latitude  as  well  as  altitude.  A  prominent  feature  common  to 
each  of  the  three  species  is  that  higher  concentrations  are 
always  found  in  the  tropical  stratosphere.  This  behavior  is 
expected  from  the  known  meteorological  patterns  through  which 
upward  transport  occurs,  preferentially  in  tropical  latitudes, 
followed  by  poleward  transport  within  the  stratosphere  itself. 

It  is  also  evident  that  the  uv  solar  radiation  which  dis¬ 
sociates  the  above  gases  is  signif icnatly  more  intense  in  the 
tropical  stratosphere  than  in  the  mid-latitude  stratosphere. 
Given  the  strong  latitudinal  variations  in  both  the  mixing 
ratios  and  the  photolysis  rates,  it  seems  natural  to  inquire 
(1)  whether  the  lifetimes  predicted  by  current  1-D  (vertical 
only)  models  are  realistic,  and  (2)  whether  the  2-D  model 
together  with  the  field  data  can  be  used  in  future  studies 
to  improve  the  lifetime  estimates? 

Using  our  2-D  models,  we  calculated  the  atmospheric  life¬ 
times  of  these  gases  and  examined  the  effect  of  the  latitudinal 
variations  of  the  mixing  ratios  and  photolysis  rates.  The 


results  have  been  published  (Ko  and  Sze,  1982).  A  brief 
discussion  will  be  given  below. 

The  atmospheric  lifetime  (t)  of  a  trace  gas  that  has 
attained  a  steady  state  concentration  in  the  atmosphere  can 
be  defined  as 


/  global  burden  \  _  /  /fNdV  \ 
\global  removal/  \/LNfdv/ 


(4.2) 


where  N  is  the  air  number  density,  f  and  L  are  the  volume 
mixing  ratio  and  loss  frequency  of  the  trace  gas  respec¬ 
tively,  and  the  bracket  represents  an  annual  average  over 
time.  Formally,  a  knowledge  of  the  spatial  and  temporal 
behavior  of  f  and  L  are  required  for  calculating  t. 

The  global  burden  of  ^0#  CFC-11  and  CFC-12  may  often 
be  estimated  from  the  observed  concentrations  at  a  few  ground 
stations  since  the  bulk  of  these  gases  resides  in  the  tropo¬ 
sphere  where  they  are  well  mixed.  The  calculation  of  the 
global  removal  rate  is,  however,  less  straight  forward  as  it 
involves  the  integral  of  the  product  of  f,  L  and  N.  The 
major  contribution  to  the  integral  arises  from  the  strato¬ 
spheric  region  where  the  quantities  f,  L  and  N  vary  signif¬ 
icantly  with  altitudes  and  latitudes.  A  common  approach  is 
to  approximate  the  integral  (which  may  be  interpreted  as  a 
spatial  and  temporal  average  of  the  product)  by  the  product 
of  the  corresponding  averaged  quantities.  This  approach 
clearly  ignores  the  correlation  between  the  quantities  f, 

L  and  N  and  hence  its  contribution  to  the  integral.  We 


92 


shall  show  in  what  follows  that  the  global  removal  rate 
depends  not  only  on  the  vertical  structure  of  f,  L  and  N  in 
the  stratosphere,  but  also  on  their  latitudinal  distributions 
and  correlations  as  well. 

The  distribution  of  f  calculated  from  the  two-dimensional 
zonal-mean  model  is  representative  of  the  zonal  average  concen¬ 
trations.  If  one  assumes  that  the  zonal  deviation  of  f  is 
small,  the  global  removal  rate  can  accurately  be  approximated 
by  integrating  the  product  of  the  zonal-mean  f  and  zonal-mean 
L. 

Figure  4-33  shows  the  calculated  column  removal  rate  of 

N20  as  a  function  of  time  and  latitude  from  our  2-D  model. 

Note  the  strong  latitudinal  gradient,  a  feature  reflecting 

the  important  effect  of  the  latitudingal  variations  in  local 

concentration  and  local  solar  flux.  During  Spring  and  Fall 

seasons,  the  maxima  occur  around  the  equator,  whereas  during 

Summer  and  Winter  seasons,  the  maxima  occur  about  15°  above 

and  below  the  equator  respectively.  Our  calculation  shows 

that  the  net  contribution  to  the  global  removal  rate  from 

the  covariance  of  f  and  L  is  positive.  The  annual  average 

9  -2-1 

column  removal  rate  for  N20  is  1.2x10  molecules  cm  s 

18  —2 

With  a  calculated  column  abundance  of  6  x  10  molecules  cm  , 
a  lifetime  of  159  years  is  derived.  This  value  may  be  com¬ 
pared  with  values  of  175  and  150  years  calculated  by  Johnston 
et  al.  (1979)  and  Levy  et  al.  (1979)  respectively. 

The  column  removal  rate  of  CFC-12  has  features  very 
similar  to  that  of  N20.  For  CFC-11,  the  latitudinal  gradient 
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Figure  4-33 

Calculated  N-O  column  removal  rate  as  a  function 
of  latitude  and  time  of  the  year. 
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of  the  removal  rate  is  somewhat  stronger  as  shown  in  Figure 
4-34.  The  calculated  lifetimes  are  65  and  135  years  for 
CFC-11  and  CFC-12  respectively.  The  corresponding  2-D  model 
calculated  lifetimes  reported  by  Hinds  (1979)  are  81  and  102 
years  and  those  from  the  Oxford  University  models  (DOE,  1979) 
are  78  and  133  years  respectively. 

The  question  of  whether  these  calculated  values  truly 
represent  the  actual  atmospheric  lifetimes  can  be  answered 
only  following  further  validation  of  the  corresponding  cal¬ 
culated  distributions  using  more  observational  data  and 
improvement  in  transport  and  chemical  data.  However,  we 
believe  that  the  latitudinal  behavior  of  the  column  removal 
rate  simulated  by  the  2-D  model  is  qualitatively  correct  and 
that  the  positive  correlation  between  the  latitudinal  varia¬ 
tion  of  f  and  L  should  not  be  ignored. 

Our  analysis  points  to  two  limitations  of  1-D  models  in 
regards  to  the!  xifetime  predictions.  First  is  their  in¬ 
ability  to  incorporate  the  covariance  effects  discussed  here. 
Second,  and  more  problematic,  is  the  inherent  difficulty  in 
validating  their  lifetime  predictions  using  field  data.  Our 
discussion  points  to  the  fact  that,  even  for  a  fully  valid¬ 
ated  1-D  model  in  which  the  calculated  profiles  are  in  agree¬ 
ment  with  global  average  observations,  the  model  may  still 
underestimate  the  global  removal  rate  and,  therefore,  over¬ 
estimate  the  lifetime  because  of  its  inability  to  account 
for  the  positive  contribution  from  the  covariance  of  f  and 
L.  With  a  modest  increase  in  computer  time,  2-D  models,  on 
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the  other  hand,  can  simulate  the  latitudinal  and  seasonal 
variations  of  atmospheric  species  as  well  as  the  solar  inso¬ 
lation  and  may  therefore  be  useful  in  refining  the  lifetime 
estimates  calculated  by  1-D  models. 

4.4  Summary  of  Model  Developments 

In  addition  to  the  features  described  in  Appendix  A, 
the  following  additional  capabilities  were  implemented  in 
our  model  during  this  past  year. 

In  the  previous  contract  year,  a  diurnal  code  was  devel¬ 
oped  for  the  1-D  model  (Sze  et  al.,  1980)  for  calculating 
the  radical  species  concentration.  With  the  grouping  tech¬ 
nique,  the  partition  of  the  radical  species  is  determined 
by  local  photochemical  reactions  independent  of  transport. 
Thus,  the  same  code  can  be  adopted  for  the  2-D  model  with 
minor  modification.  The  2-D  model  has  the  capability  of 
calculating  the  diurnal  behavior  of  all  the  short-lived 
species  once  the  concentrations  of  the  long-lived  species 
are  specified.  The  results  of  the  calculations  are  being 
prepared  for  publication.  In  addition,  in  the  past  month, 
we  have  explroed  ways  to  include  diurnal  effects  in  the  03 
equation. 

The  present  2-D  model  uses  an  explicit  time  scheme  to 
solve  the  trace  gas  equation.  To  ensure  numerical  stability, 
a  4-hour  time  step  has  to  be  used  to  propagate  the  solutions 
in  time.  We  have  succeeded  in  devising  a  semi-implicit  scheme 
using  the  alternative  direction  method.  With  the  new  scheme. 


the  time  step  can  be  taken  to  be  1  day.  Further  work,  how¬ 
ever,  is  required  to  apply  the  method  to  the  full  model. 

Finally,  in  the  past  year,  we  have  been  assessing  alter¬ 
nate  methods  of  treating  the  dynamic  transports  in  the  model 
so  as  to  include  the  coupling  between  radiation-dynamics  and 
chemistry.  The  Eulerian  formalism  used  in  the  present  models 
required  values  for  eddy  diffusion  coefficients  which  are 
parameterized  from  observations  of  the  present  day  atmosphere. 
Among  the  recently  suggested  alternate  methods,  we  found  the 
formalism  of  residual  circulation  to  be  most  attractive 
(Holten,  1981;  Tung,  1982) .  We  are  at  present  exploring 
ways  to  test  some  of  these  alternative  approaches. 

Although  our  2-D  model  is  performing  to  our  expectations, 
we  still  feel  some  of  these  improvements  will  add  to  the  capa¬ 
bility  of  the  model.  We  will  be  looking  for  alternate  sources 
of  funding  for  this  work. 


98 


5.  EFFECT  OF  TRACE  METAL  ON  STRATOSPHERIC  CHEMISTRY 


5 . 1  Background 

It  is  now  believed  that  trace  metals  are  introduced 
into  the  earth's  atmosphere  via  ablation  of  meteors.  The 
metal  diffuses  downward  and  is  found  in  neutral  form  around 
and  below  90  km.  Meteors  contain  various  metals  including 
about  10%  by  weight  each  of  Fe  and  Mg,  with  Na,  Ca,  Al  and 
Ni  each  constituting  about  1-2%  by  weight.  Among  the  metals, 
Na  is  most  studied  because  of  the  observation  of  the  Na  layer 
around  90  km  as  deduced  from  spectroscopic  (Cf.  see  Hun ten, 
1967)  and  laser  radar  measurements  (Sanford  and  Gibson, 

1970;  Blamont  et  al.,  1972;  Megie  and  Blamont,  1977,  Simonich 
et  al. ,  1979) . 

The  basic  sequence  of  reactions  responsible  for  main¬ 
taining  the  neutral  sodium  concentration  was  first  suggested 
by  Chapman  (1939)  to  be 


*1 

Na  +  03 

NaO  +  02 

NaO  +  0 

* 2 

Na(2P)  +  0 

However,  since  there  is  no  laboratory  measurement  of  the 
kinetic  reaction  rates,  it  can  only  be  demonstrated  that 
with  reasonable  assumed  rates  for  k^  and  k2  and  observed  0 
and  03  concentrations,  one  can  deduce  distributions  for  Na 
that  are  consistent  with  observations. 


99 


Recently,  Ferguson  (1978)  argued  that  the  mass  spectro- 
metric  measurement  of  positive  ions  in  the  stratosphere 
(Arnold  et  al.,  1977,  1978)  can  be  interpreted  as  NaOH 
clusters.  Furthermore,  he  estimated  the  NaOH  concentrations 
to  be  ^103  cm  3  at  40  km.  The  modeling  work  of  Liu  and  Reid 
(1979)  was  one  of  the  more  recent  studies  that  extended  the 
calculation  of  the  sodium  species  down  to  the  stratosphere. 
Using  the  reaction  scheme  summarized  in  Table  5-1,  they  suc¬ 
ceeded  in  simulating  the  neutral  sodium  layer  at  90  km  with 
a  concentration  of  3  x 10  cm  J  and  found  NaOH  to  be  the  most 

abundant  sodium  species  below  80  km.  In  addition,  the  NaOH 

.  5  -3 

concentration  at  40  km  is  calculated  to  be  ^3 x  10  cm  ,  in 

agreement  with  Ferguson's  estimates. 

Murad  and  Swider  (1979)  re-examined  the  thermodynamic 

data  on  the  sodium  species  and  pointed  out  the  following 

reactions  are  estimated  to  be  exothermic  and  therefore  might 

have  to  be  considered  in  the  sodium  chemistry  scheme. 


In  addition. 


Na02  +  0 

-► 

NaO  +  02 

k13 

NaO  +  03 

-► 

Na  +  20 2 

k14 

NaO  +  H20 

NaOH  +  OH 

k15 

NaOH  +  Cl 

-► 

NaCl  +  OH 

k16 

NaOH  +  CIO 

NaCl  +  H20 

k17 

NaOH  +  HC1 

-*• 

NaCl  +  H20 

k18 

the  reaction  of  NaOH  with  HN03  and  C02, 
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Reactions 


Rate  Constants 


(a) 


K1 

Na  +  03 

Na  0  4  02 

3.3 

X 

10"10 

(b) 

K2 

NaO  +  0 

-* 

Na  +  02 

1.6 

X 

10"10 

(b) 

K3 

Na  +  02  +  M 

-► 

Na02  +  M 

8.4 

X 

IQ’34 

exp(2|0)  (c) 

K4 

NaO  +  03 

Na02  +  02 

5 

X 

ID'11 

K5 

NaO  +  H02 

->■ 

NaOH  +  02 

1 

X 

ic-11 

K6 

Na02  +  H 

-► 

NaOH  +  0 

7 

X 

lO"12 

K7 

Na02  +  OH 

NaOH  +  02 

1 

X 

io‘n 

KSa 

Na  +  H02 

■+ 

NaOH  +  0 

7 

X 

10~12 

K8b 

Na  +  H02 

NaH  +  02 

1 

X 

10'11 

K8c 

Na  +  H02 

NaO  +  OH 

8 

X 

10-12 

k9 

NaOH  +  hv 

-► 

Na  +  OH 

same  as  photolysis 
of  water 

K10 

NaOH  +  0(1D) 

-► 

NaO  +  OH 

2.3 

X 

10-10 

K12 

NaOH  +  H 

-► 

Na  +  H20 

2 

X 

10"11 

.-800, 
exp  (  -^~) 

3—1  6  — 

(a)  Units  are  cm  s  for  two-body  reactions  and  cm  s 

for  three-body  reactions.  All  rates  are  assumed 
except  as  noted  below. 

(b)  Estimated  by  Kolb  and  Elgin  (1976) . 

(c)  Estimated  using  high  T  measurement  of  Carabetta  and 
Kashan  (1968) . 
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-►  NaN03  +  H20 


NaOH  +  HN03 
NaOH  +  C02  +  M  -*•  NaHC03  +  M 


could  conceivably  be  important,  although  there  is  a  lack  of 
data  for  NaN03  and  NaHC03  to  estimate  the  energetics  for  the 
reactions. 

Of  the  reactions  suggested,  reactions  16  through  20  are 
of  particular  interest  if  NaOH  is  the  dominant  sodium  species 
in  the  stratosphere.  In  a  subsequent  paper,  Murad  et  al. 
(1981)  argued  that  other  metals  from  meteors  may  proceed  via 
analogous  reactions  to  combine  with  chlorine  species  to 
form  metal  chlorides  which  may  thon  be  incorporated  in  aero¬ 
sols  arid  removed  via  depositions.  It  was  further  suggested 
that  high  flying  aircrafts  may  be  another  source  of  trace 
metals  through  combustion  of  jet  fuel  and  engine  erosion. 

Our  study  reported  here  was  initiated  to  assess  the  impact 
of  the  reaction  of  the  trace  metals  as  a  sink  for  chlorine 
species  in  the  stratosphere  with  particular  attention  to  the 
possible  influence  on  the  chlorine-catalyzed  destruction  of 


5.2  Model  Study 

We  first  used  our  1-D  model  to  estimate  the  steady  state 
concentrations  of  total  metal  in  the  atmosphere.  The  trace 
metals  are  either  injected  at  the  mesopause  to  simulate  down¬ 
ward  diffusion  of  meteoric  metal  from  the  mesosphere  or  in¬ 
jected  in  the  lower  stratosphere  to  simulate  aircraft  emission. 
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In  all  cases,  a  tropospheric  sink  with  a  lifetime  of  10  day 
is  used  to  simulate  rainout  and  deposition. 

Figure  5-1  shows  the  calculated  total  metal  concentra¬ 
tions  as  a  function  of  altitude  for  several  cases  with  dif¬ 
ferent  injection  parameters: 

6  -o 

(i)  top  injection  of  2x10  metal  molecules  cm 

-1  4 

s  (2x10  ton/yr)  at  71  km  to  simulate  the 
effect  of  ablation  from  meteor  (curve  marked  Mx) 

3 

(ii)  injection  of  2  x  10  ton/yr  of  metal  spread 
uniformly  over  a  1  Jem  thick  band  at  16  Jon  and 
20  Jem  respectively  to  simulate  effect  of  air¬ 
craft  engine  emissions  for  different  modes  of 
operation  (curve  marked  I1(j  and  I2Q) 

Studies  (CIAP ,  1975a, b)  indicated  that  the  emissions 
of  metals  from  engine  erosion  is  much  higher  than  those 
associated  with  fuel  combustion.  The  metal  content  in  jet 
fuel  is  estimated  to  be  in  the  ppb  level  while  metals  com¬ 
ponent  of  exhaust  from  engine  erosion  could  be  in  the  ppm 

.  9 

level.  Using  the  CIAP  reported  estimate  of  3x10  ton/yr 

for  the  total  exhaust  associated  with  aircraft  operations, 

we  arrived  at  the  value  of  2  x  10"5  ton/yr  by  assuming  that 

the  exhaust  contains  ^1  ppm  of  metal. 

6  —  2  “1 

The  value  of  2  x 10  cm  s  for  meteoric  input  is 

4 

based  on  Gadsden's  (1976)  estimate  of  9x10  ton/yr  for 

total  meteor  flux  and  that  about  20%  of  the  meteor  flux  is 

ablated  to  produce  metal  atoms.  The  estimates  for  total 

2  7 

meteoric  flux  range  from  10  to  10  tons/yr  (Hughes,  1972) 


MIXING  RATIO  (ppbv) 


Figure  5-1 

Calculated  total  metal  profiles  for  different  injection 
rates.  See  text  for  injection  parameters.  The  C1X 
profile  is  included  for  comparison. 


4 

although  values  of  order  10  ton/yr  are  most  often  used. 

6  -2  —I 

It  should  be  noted  that  the  value  of  2  x  10  cm  s  is 

2 

about  a  factor  of  10  larger  than  the  values  of  Na-flux 
used  by  Liu  and  Reid  (1979) .  Given  the  fact  that  about  a 
few  percent  of  the  metal  flux  is  in  the  form  of  Na,  our 
assumed  input  rate  seems  reasonable.  We  will  now  restrict 
our  discussion  to  meteor  input  case  since  the  total  metal 
concentration  resulting  from  aircraft  emissions  are  much 
smaller. 

In  this  study,  following  Murad  et  al.  (1981) ,  we 
assumed  that  once  the  metal  chloride  (MCI)  is  formed,  it 
will  be  removed  via  heterogeneous  processes  and  that  there 
is  no  pathway  to  recycle  MCI  back  to  free  metal  and  chlorine. 
The  calculated  total  chlorine  (C1X)  in  the  present  day  atmo¬ 
sphere  is  also  plotted  in  Figure  5-1.  In  the  case  of  no 
recycling,  each  metal  atom  can  only  react  with  the  chlorine 
species  once.  Since  there  is  a  -comparable  amount  of  C1X  and 
MX  above  50  km,  most  of  the  C1X  can  be  converted  to  MCI  if 
the  reaction  is  sufficiently  fast.  In  contrast,  the  MX  con¬ 
centration  is  about  1/3  of  the  C1X  concentration  around  40 
km.  The  reaction  between  metal  and  chlorine  species  can  at 
most  reduce  the  free  chlorine  concentration  by  1/3  even  if 
conversion  is  complete. 

The  amount  of  free  metal  (MY  =  MX  -  MCI)  in  the  atmo¬ 
sphere  will  clearly  depend  on  whether  the  rate  of  conversion 
is  slow  enough  for  the  free  metal  to  diffuse  downward.  Since 
HC1  is  the  most  abundant  chlorine  species,  the  conversion 
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rate  for  MX  to  MCI  can  be  estimated  using  the  reaction 

MOH  +  HC1  -*■  MCI  +  H20  k21  . 

Figure  5-2  shows  the  mixing  ratio  for  MY  for  different 

assumed  rates  of  k21.  Murad  et  al.  (1981)  estimated  that 

-9  3-1 

reaction  21  can  proceed  as  fast  as  10  cm  s  Thus, 
there  should  be  very  little  free  metal  in  the  atmosphere. 

Figure  5-3  shows  the  distribution  of  the  chlorine 
species  for  the  case  with  k21 = 10  11cm3  s_1.  Note  that 
although  the  free  chlorine  concentration  is  reduced  signif 
icantly  above  50  km,  there  is  little  impact  below  40  km 
because  of  lack  of  metal  atoms.  The  chlorine  catalytic 
destruction  of  ozone  is  most  effective  around  35  km.  Thus 
the  overall  impact  of  metal  on  the  chlorine-ozone  cycle  is 
small  unless  we  have  grossly  underestimated  the  amount  of 
metal  in  the  stratosphere. 

5.3  Recent  Developments 

Since  our  study,  there  have  been  several  new  reports 
on  the  kinetic  data  on  several  of  the  reactions.  Recent 
work  by  Husain  and  Plain  (1982)  indicates  a  rate  of  10 
cm6  s  1  for  k^-  Most  recently,  Rowland  and  Makide  (1982) 
found  NaOH  absorbs  strongly  in  the  ultraviolet  region  and 

3 

has  a  photolysis  rate  of  about  a  factor  of  10  faster  than 
those  assumed  in  Table  5-1.  With  these  new  data,  Na02 
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Figure  5-2 

Calculated  mixing  ratio  profiles  for  free  metal 
(MY  =  MX  -  MCI)  for  different  assumed  values  of 
reaction  rate  k-,. 


Case  A 
Case  B 
Case  C 
Case  D 


1<T  ^  cml  s~j 
10  „  cm,  s't 
10  1  cmJ  s"1 
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became  the  dominant  species  below  80  km  with  NaOH  playing 
a  secondary  role. 

With  regard  to  the  positive  ion  clusters  data,  Arnold 
and  Henschen  (1981)  reported  that  they  failed  to  detect  any 
Na+  thus  negating  Ferguson's  (1978)  estimates  of  NaOH  con¬ 
centration. 

Another  recent  development  concerns  the  atmospheric 
fate  of  NaCl.  Recent  work  by  Rowland  and  Rogers  (1982) 
indicates  that  NaCl  photolysizes  rapidly  and  has  strato¬ 
spheric  lifetimes  of  about  a  few  minutes.  Reaction  k^g 
followed  by  photolysis 

NaCl  +  hv  -*■  Na  +  Cl  k»„ 


is  equivalent  to 


NaOH  +  HC1  •*  Na  +  H20  +  Cl 

and  may  lead  to  a  potential  catalytic  pathway  for  converting 
HCl  to  Cl.  Once  the  catalytic  cycle  is  introduced,  the  im¬ 
pact  of  metal  is  no  longer  limited  by  its  concentrations. 

It  seems  prudent  to  re-examine  the  role  of  metal  chemistry 
in  regard  to  other  chemical  cycles  in  view  of  the  new 
reports  on  kinetic  data. 
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A.l  Introduction 


The  AER  zonal-mean  model  is  designed  to  simulate  the 
zonal-mean  behavior  of  atmospheric  trace  gases.  The  domain 
of  the  model  covers  the  globe  from  pole  to  pole  and  up  to 

an  altitude  of  60km.  The  horizontal  and  vertical  resolu- 

.  .  3 

tions  are  approximately  10  km  and  3.5km  respectively.  The 
dynamic  variables  are  governed  by  the  system  of  equations 
obtained  by  taking  the  Eulerian  zonal-mean  of  the  simplified 
primitive  equations  which  assumes  geotrophic  wind  balance. 
Dynamic  transport  of  trace  gases  is  affected  by  advection 
from  zonal-mean  wind  and  eddy  transport  parameterized  by  a 
symmetric  diffusion  tensor.  The  photochemical  scheme  used 
in  the  2-D  model  is  based  on  the  grouping  techniques  and 
calculates  over  35  species  affecting  the  03  balance.  A 
24-hour  averaging  of  the  photolysis  rate  calculation  is  used 
to  simulate  diurnal  average  condition.  The  behavior  of  the 
trace  gases  in  the  atmosphere  is  governed  by  the  trace  gas 
equation 


N  (- 


3£i 

3t 


+  v 


vfi) 


w 


(A.l) 


where  N  is  the  bulk  air  density,  v  is  the  velocity  field  of 
bulk  motion,  f^  is  the  volume  mixing  ratio  and  the  local 
net  production/loss  rate  of  the  ifc^  species. 

Equation  (A.l)  is  to  be  considered  as  a  system  of  equa¬ 
tions  giving  the  time  rate  of  change  of  the  Eulerian  field 
variables  f^'s  as  functions  of  the  spatial  coordinates  in 


in  three  dimensions.  To  solve  the  system,  one  requires 
knowledge  of  N,  which  depends  on  atmospheric  conditions, 
pressures  and  temperature:  the  velocity  wind  field  v  and 
knowledge  of  the  chemical  interactions.  To  reduce  the 
equations  to  two  dimensions,  a  method  for  averaging  must 
be  specified.  Each  of  these  will  be  discussed  in  the  follow' 
ing  sections. 

The  structure  of  the  dynamic  framework  in  the  model 
follows  closely  that  of  the  Oxford  model  (Harwood  and  Pyle 
1975) .  The  Oxford  model  was  originally  designed  to  calcu¬ 
late  the  dynamic  variables  together  with  the  trace  gas 
concentrations.  At  the  present  stage  of  the  development, 
the  AER  model  uses  pre-specified  dynamic  variables  as 
input  and  solve  only  the  trace  gas  equations.  The  assump¬ 
tions  used  in  arriving  at  the  system  of  dynamic  equations 
in  the  zonal  model  will  be  reviewed  in  Section  A. 2  and  A. 3. 
The  discussion  on  the  treatment  and  method  of  solution  of 
the  system  of  dynamic  equations  is  kept  to  a  minimum.  The 
reader  is  referred  to  in  the  original  paper  of  Harwood  and 
Pyle  (1975) .  The  treatment  of  the  trace  gas  equations  will 
be  discussed  in  more  detail  in  Section  A. 4  through  Section 
A. 7.  The  overall  model  algorithm  will  be  discussed  in 
Section  A. 8. 

A. 2  The  System  of  Equations 

The  atmospheric  circulation  is  governed  by  the  coupled 
system  of  dynamic  and  thermodynamic  equations  (cf.  Lorentz, 


124 


1967)  : 


momentum  equation 

H  -  -22*y  -  m  »  -  J* 


(A. 2) 


thermodynamic  equation 


d6 

dt 


J 


(A. 3) 


continuity  equation 


dN 

dt 


-NV- v 


(A. 4) 


equation  of  state:  ideal  gas  law 


P  *  NRT  (A.  5) 

d  3 

where  ^  ^  +  y  •  V  is  the  total  time  derivative 

ft  =  angular  velocity  of  Earth 

M  =  average  mass  of  an  air  molecule 

p  =  pressure 

4>  *  the  geopotential  gz  where  g  is  acceleration  due  to 
gravity,  z  is  geometric  altitude 

6  =  potential  temperature  related  to  temperature 

T  by  6  =  T  ^  with  K  =  R/Cp?  R  the  gas  constant 
and  Cp  the  specific  heat  at  constant  pressure 

T 

J  =  the  diabatic  influence  with  C  J*-  the  heating  rate 
per  unit  mass  p 


R  =  the  gas  constant 

T  =  temperature 

Note  that  frictional  forces  have  been  left  out  in  the  momen¬ 
tum  equation  (A. 2)  under  the  assumption  that  they  are  unimpor¬ 
tant  for  large  scale  motion. 

The  system  of  equations  (A. 2)  to  (A. 5)  is  coupled  to  the 
trace  gas  equation  through  the  J  term  which  depends  on  dis¬ 
tribution  of  gases  such  as  0^,  C02 ,  N20  and  CH^  in  the 
atmosphere.  Thus,  in  principle,  equations  (A.l)  through  (A. 5) 
should  be  solved  simultaneously  as  a  system.  The  following 
physical  assumptions  are  usually  adopted  to  obtain  the  so- 
called  primitive  equations  (cf .  Holten,  1375) : 

A)  Replacement  of  the  vertical  momentum  equation  by 
the  hydrostatic  equilibrium  condition 

§f  -  -MN  g  (A. 6) 

B)  Discard  terms  in  the  expression  for  the  coriolis 
force  that  contain  the  vertical  component  of  the 
velocity. 

C)  Replace  r  in  the  resulting  equation  by  the  mean 
distance  from  the  earth's  center. 

In  order  to  write  the  momentum  equation  in  scaler  form,  it 
is  convenient  to  introduce  a  coordinate  system  where  the 
pressure  p  is  used  as  the  vertical  coordinates  in  conjunction 
with  X  the  longitude  and  <p  the  latitude.  In  this  coordinate 
system,  the  components  of  the  velocity  vectors  are; 
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(A.  7) 


=  a  cos< 


With  the  hydrostatic  equilibrium  assumption,  the  continuity 
equation  and  the  trace  gas  equation  take  the  form 


1  3u  ,  1  d  ,  . .  3w 

r5 55?  5T  +  —cob*  5*  (v  «>■♦)  +  jp  -  « 


(A.  8) 


oi.  13  Q 

TT  +  rssi?  5X  <fiu»  +  T-kz*  n  <V  «»♦>  +  -k  ,fiw>  - 

(A. 9) 

The  primitive  equations  can  be  written  (cf.  Lorentz, 
1967;  Holten ,  1975) 

3u  .  1  3,2..  1  3  a 

3t  +  a  cos<fr  3X  (u  }  +  a~co8(|>  3?  (uv  cos*)  +  ab  (uw> 


g-v  *an*  -  2(lv  sin, - 1—  |*  .  o 

a  a  cos4>  3X 


9t  a  cos 4>  3X  +  a  cos4>  3 4>  cos$)  +  jjp  (vw)  (A. 11) 
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(A. 12) 
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(A. 13) 
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(A. 14) 
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(A.  15) 


2 

In  the  above  equations,  the  terms  and  - — are 

A  A 

curvature  terms  arising  from  the  non-Euclidean  nature  of 

1  84 

the  coordinate  frames.  The  geopotential  term  —  and 
1  84 

-  s  i  -gy  has  been  written  out  explicitly  with  4  ■  g*. 


In  this  arrangement,  equations  (A. 9)  through  (A. 12) 
are  prognostic  equations  for  f,  u,  v  and  6  respectively, 
while  equations  (A. 13)  and  (A.  14)  can  be  considered  as 
diagnostic  equations  for  w  and  N.  Since  p  is  being  used 
as  an  independent  variable,  equation  (A. 15)  serves  the  pur¬ 
pose  of  transforming  between  geometrical  altitude  and  the 
pressure  coordinate. 


A.  3  Zonal  Averaging 

We  will  now  discuss  ways  of  reducing  the  system  of 
primitive  equations  to  2  dimensions  for  numerical  modeling. 
Given  any  of  the  dependent  Eulerian  variables  in  the  primi¬ 
tive  equations,  e.g.,  f,  by  averaging  over  one  of  the  spatial 
coordinates,  one  obtains  f  (t,y)  where  y  is  the  vector 
representing  the  two  remaining  coordinates. 
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The  corresponding  equations  governing  the  averaged 
variables  can  be  obtained  by  applying  the  averaging  opera¬ 
tion  to  the  primitive  equations.  In  atmospheric  modeling 
studies,  one  is  interested  in  averaging  over  the  longitud¬ 
inal  coordinate  to  obtain  £  as  a  function  of  time,  latitude 
and  altitude,  which  can  be  viewed  as  the  zonal-mean  concen¬ 
tration.  However,  the  exact  physical  interpretation  of  f 
must  depend  on  the  averaging  process  involved. 

In  the  Eulerian-mean  model,  the  averaging  process  is 
performed  over  Eulerian  field  quantities  by  integrating 
over  the  longitude  while  holding  all  other  coordinates 
fixed.  Thus,  given  g(t,4>,p,X),  (where  t  is  time,  $  the 
latitude,  p  the  pressure  height  coordinates  and  X  the  longi 
tude) ,  one  obtains 

i(t,o,P)  =  f- (A. 16) 

It  is  convenient  to  define  g',  the  deviations  from  the 
mean ,  by 

g'  (t,$,p,X)  «  g  (t,<J>,p)  -  g(t,4>,p,X)  (A. 17) 

It  follows  from  equation  (A. 16)  that 


g' (t,$,p,A)  =  0 


(A. 18) 


and 


gK  -  gh  +  g ' h ' 


(A. 19) 
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Before  applying  the  Eulerian-mean  operation  to  the 
primitive  equations,  the  geotrophic  wind  balance  assumption 
will  be  used  to  replace  (A. 11)  by 

2flusin<J>  "  |  -  0  (A. 20) 


Using  the  fact  that  the  " - H  operation  commutes  with  the 

coordinate  differential  operators,  we  have 
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(A. 22) 
(A. 23) 
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(A. 27) 
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Equations  (A. 21,  A. 22  and  A. 24)  are  to  be  considered  as 
prognostic  equations  for  f,  u  and  0  respectively.  Equation 
(A. 25)  implies  existence  of  a  stream  function  ij;,  where 


v  = _ L_  11 

COS<f>  dp 

(A. 28) 

W  =  -  1  -  11 

a  cos  <j>  3  <f> 


By  differentiating  w.r.t. 
equation  (A. 2 7),  Equation 
thermal  wind  relation 


3u 

3p 


(J)  (S) 

2flsin4> 


1 

a 


p  and  using  the  hydrostatic 
(A. 2 3)  can  be  rewritten  as  the 


36 


(A. 29) 


To  obtain  an  equation  for  ty,  one  notes  that  (A. 22)  and  (A. 24) 
is  of  the  form 


3u 

3t 
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(A. 30) 


il 

3t 


B 


(A.  31) 


where  A  and  B  contains  only  spatial  derivatives.  By  differ¬ 
entiating  (A. 30)  w.r.t.  p  and  (A. 31)  w.r.t.  $  and  using  the 
thermal  wind  relation,  we  obtain 
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Conversion  to  geometrical  altitude; 


3z  RT 

Tfp  ~  ”  Mgp 


(A. 39) 


Given  the  variables  f.,  u,  T,  at  time  t,  if  the  sources 
Qi 

terms  (-jg-)  .•  Ff .  #  J#  FQ  and  Fu  are  known,  (A. 36)  can1  be  used 
to  solve  from  ip.  The  velocity  field  v  and  w  are  specified 
by  (A. 37).  Equations  (A. 33) ,  (A. 34)  and  (A. 35)  can  then  be 
used  to  solve  for  f^  u  and  T  at  a  later  time  t  +  At.  Thus 
the  whole  system  can  be  propagated  forward  in  time. 

A. 4  Source  Terms  for  Zonal-Mean  Trace-Gas  Equations 

We  shall  now  restrict  our  discussion  to  the  treatment 

of  the  trace-gas  equations.  A  discussion  on  the  treatment 

of  the  dynamic  equations  can  be  found  in  Harwood  and  Pyle 

(1975).  Solution  to  (A. 33)  depends  on  knowledge  of  v,  w, 

Qi 

(-jj-)  and  Ff^.  As  pointed  out  in  the  introduction,  v  and 
v  are  obtained  from  input  files.  The  term  Ff^  is  parameter¬ 
ized  by  mean  of  an  eddy  diffusion  tensor  (Reed  and  German, 
1965) .  In  the  present  coordinate  system, 


(A. 40) 


where  =  K^.  Note  that  it  is  assumed  that  the  same 
diffusion  tensor  can  be  used  for  all  species. 
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The  formulation  of  Reed  and  German  (1965)  is  based  on 
mixing  length  type  argument.  Let  the  motion  in  a  fluid  be 
given  by  a  mean  motion  v  with  deviation  v* .  For  a  quantity 
x  which  is  conserved  along  the  mean  flow,  it  is  reasonable 
to  argue  that  the  deviation  from  the  mean  x'  is  given  by 

x'  =  -  1'  •  Vx 

where  1*  is  the  displacement  vector  associated  with  the 
deviation  velocity  v' .  If  one  further  assumes  that  1*  and 
v'  are  linearly  related,  i.e. 

1'  =  v'At 

where  At  is  the  characteristic  time-scale  of  the  deviation 
from  the  mean  motion,  it  follows  that 


v'x'  *  -  v*  U'-  Vx)-  -  At  l'l '  •  Vx 


(A. 41) 


where  XTTT  i*  a  symmetry  tensor.  In  scaler  form,  we  have 


(A. 42) 


If  one  assumes  that  the  eddy  fluxes  of  potential  temperature 
can  be  parameterized  this  way,  one  can  deduce  the  values  for 
the  tensor  component  by  examining  the  temperature  and  wind 
field  statistics  (cf.  Reed  and  German,  1965;  Luther  1973). 
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The  exact  form  of  the  term  will  be  discussed  in 
Section  A. 6.  The  discussion  here  will  be  restricted  to  the 
treatment  of  the  averaging  processes.  The  term  Q/N  consists 
of  sums  of  terms  each  of  which  take  one  of  the  following 
forms 


Ji£i!  Nkij(T>  fi V  Nkijk(T)£ifi£fc 


for  photolytic  processes  and  two-body  and  three-body  reac¬ 
tions  respectively.  Note  that  the  temperature  dependence 
of  the  reaction  rates  k • .  and  k. ..  are  explicitly  displayed. 

i  j  1 

The  longitudinal  behavior  of  J.  and  f  arise  from  the  diurnal 

1  i 

effect  since  each  longitude  is  at  a  different  local  time 
at  any  instant  in  time.  In  the  present  model,  the  term 
(®)  is  formed  by  taking  the  corresponding  products  of  the 
zonal-mean  quantities.  No  attempt  is  made  to  take  into 
account  for  the  covariance  of  the  individual  terms  making 
up  the  product. 


A.  5  Numerical  Solution  of  the  Trace-Gas  Equations 

In  finding  the  numerical  solution  of  the  equation,  it 
is  more  convenient  to  use  the  log-pressure  coordinate 
£  =  ln(pQ/p)  where  pQ  is  taken  to  be  1000  mb.  In  this 
coordinate  system,  the  vertical  velocity  is  given  by 


u 


=  as. 

dt 


Equation  (A. 33)  can  be  rewritten  in  the  form: 


(A. 43) 
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K« 


kK«  K« 


K*«  -  K«  -  -  K*  e^ 

K*e  .  KPP(|1)2  „  KPP  e2^ 


To  obtain  a  numerical  solution  of  the  equations,  the 
differential  operators  are  first  converted  to  finite  differ¬ 
ence  form.  The  two-dimension  space  is  divided  into  a 
19  x  17  grid.  The  globe  is  covered  by  19  latitude  belts 

3 

each  9.5°  (10  km)  wide.  The  vertical  resolution  is  con¬ 
stant  in  £  *  In  pQ/p  with  *  0.5  or  &>z  *  3.5  km.  The  top 
boundary  corresponds  to  0.3  mb  level  or  ^60  km. 

The  grid  is  schematically  represented  in  Figure  A-l. 
Within  each  box,  the  variables  are  defined  at  different 
positions  so  that  all  the  individual  terms  entering  into 
Equation  (A. 44)  are  defined  at  the  center  of  the  grid. 

This  is  summarized  in  Table  A-l. 
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Schematic  diagram  showing  the  grid  sizes  in  the  2-D  model 
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Equation  (A. 33)  is  used  to  define  at  each  time  step. 
The  R.H.S.  is  evaluated  using  v,  w  and  K's  from  input  files, 
(^)  calculated  from  the  chemistry  scheme  and  values  of  f^'s 
at  the  current  time  step.  The  Adam-Bashford  scheme  is 


used  to  calculate  the  f^  at  the  next  time  step,  i.e.: 


f±(t  +  At) 


(1.5  (t)  -  0.5  (t  -  At)  )At 


Adjustments  are  made  at  each  time  step  to  simulate  the 
appropriate  boundary  condition  at  the  top  and  bottom  of  the 
atmosphere.  These  will  be  discussed  in  Section  A. 7. 

The  Adam-Bashford  scheme  is  a  second  order  accuracy 
in  time.  To  ensure  stability.  At  must  be  sufficiently 
small,  with  the  present  grid  size.  At  is  taken  to  be  4  hours. 


A.  6  Dynamical  Variables  Input 

The  dynamical  variables  input  necessary  for  solving 
(A. 44)  are  the  temperature  field,  the  stream  function  for 
generating  v  and  u  and  the  eddy  diffusion  coeffidients.  The 
values  for  the  components  of  the  eddy  diffusion  coefficients 
are  from  Luther  (1973) .  The  temperature  and  ^-function  are 
from  previous  calculations  of  Harwood  and  Pyle  (1977) . 


a)  Eddy  diffusion  tensor: 

The  values  of  the  symmetric  tensor  are  compiled  by 
\  Luther  (1973)  using  the  procedure  similar  to  the  one  described 
by  Reed  and  German  (1965).  The  heat  transfer,  temperatures 
and  wind  variance  data  of  Oort  and  Rasmussen  (1971)  are 
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Figure  A- 2 

Latitude-altitude  crose-eection  plots  of  K  for  different 
seasons.  Values  are  from  Luther  (1973) . 
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Figure  A- 3 

Latitude-altitude  crosa-section  plots  of  K  for  different 
seasons.  Values  are  from  Luther  (1973) .  “ 
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Figure  A-4 


Latitude-altitude  cross-section  plots  of  K  for  different 
seasons.  Values  are  from  Luther  (1973) .  yz 


analyzed  to  derive  a  set  of  values  as  functions  of  time  and 
space.  The  original  data  only  covers  part  of  the  northern 
hemisphere.  The  results  of  Newell  et  al.  (1966)  and  Wofsy 
and  McElroy  (1973)  are  used  to  extrapolate  the  values  at 
places  where  there  is  no  data.  The  values  for  K^,  K*^ 
and  for  the  northern  hemisphere  written  on  files  as 
function  of  latitude  and  altitude  with  a  set  of  values  for 
each  month.  The  values  of  the  southern  hemisphere  is 
obtained  by  using  the  northern  hemisphere  values  shifted 
by  six  months.  The  values  for  the  K’s  are  read  from  input 
file  and  kept  constant  over  the  whole  month.  Values  for 
K's  are  plotted  in  Figure  A-2  through  A-4  for  different 
seasons. 

b)  Temperature  and  t|>-field: 

The  values  of  T  and  4>  are  tabulated  as  functions  of 
latitude  and  altitude  around  the  year  with  one  set  of  values 
every  10  days.  A  simple  linear  interpolation  scheme  is 
used  to  obtain  the  values  for  each  time-step  as  required. 
Plots  of  the  values  of  the  variables  for  different  seasons 
are  given  in  Figure  A-5,  A-6  and  A-7. 

The  values  are  from  model  calculations  of  the  Oxford 
model  in  which  the  dynamic  equations  are  solved  together 
with  the  03  equation.  A  detail  of  the  model  calculation 
can  be  found  in  Harwood  and  Pyle  (1977) ;  only  a  brief 
summary  will  be  given  here. 
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Figure  A- 5 

Latitude-altitude  cross-section  plots  of  the  temperature  for  different  seasons 
Values  are  from  Harwood  and  Pyle  (1977). 
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(i)  Treatment  of  Ozone  equation  in  Harwood  and  Pyle  (1977) 
It  is  assumed  that  deviations  from  the  photochemical 

equilibrium  concentration  of  03 ( [03 )e) decay  with  a  chemical 
lifetime  T.  In  terms  of  the  parameters  used  in  (A. 44) 

[0-]e-[03] 
q_  =  — 5 - £_ 

°3 

In  the  1977  model  run,  [03]e  and  t  are  precalculated  and 
tabulated.  A  chemical  scheme  involving  OH  and  NOX  is  used 
where  [03)e  and  r  are  calculated  from  given  latitudinal 
independent  altitude  profile  of  OH,  NOX  and  reaction  rates. 

The  eddy  fluxes  for  03  are  parameterized  by  eddy 
diffusion  tensor.  A  zero  mixing  ratio  boundary  condition 
is  imposed  at  the  lower  boundary  to  simulate  rapid  destruc¬ 
tion  at  the  ground. 

(ii)  Treatment  of  temperature  equation  in  Harwood  and  Pyle 
(1977)  : 

The  heating/cooling  rate  is  treated  as  follows: 

Solar  heating:  Heating  by  03  absorption  using  the 
parameterization  of  Kennedy  (1964)  and  model 
calculated  03 

Thermal  radiation: 

for  troposphere  (up  to  12km) ,  constant  lapse  rate 
from  climatology  is  assumed 

1.5  K/day  0-8  km 
.75  K/day  9-12  km 

for  lower  stratosphere  (12km-24km) ,  thermal  cooling 
is  set  equal  to  solar  heating  so  that  there  is  no 
net  heating 

for  stratosphere  (above  24km)  long  wave  cooling  by 
C<>2  and  O3  each  parameterized  by  two  altitude 
dependent  coefficients. 
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Eddy  fluxes  of  potential  temperature  are  parameterized 
by  eddy  diffusion  tensor.  The  heat  flux  across  the  lower 
boundary  is  calculated  as  functions  of  wind  and  temperature 
via  the  method  of  Clarke  (1970).  The  underlying  surface  is 
assumed  to  be  at  fixed  temperature  with  infinite  heat 
capacity. 

(iii)  Treatment  of  angular  momentum  equation: 

The  momentum  flux  across  the  boundary  simulating  fric¬ 
tional  drag  is  calculated  from  the  wind  and  temperature  via 
the  method  of  Clarke  (1970).  The  values  for  horizontal 
eddy  momentum  fluxes  are  deduced  from  the  observation  of 
Nimbus-V  SCR  data.  Vertical  eddy  momentum  is  ignored. 

A. 7  Chemistry  Scheme 

The  design  of  the  chemistry  scheme  is  based  on  the 
grouping  technique.  The  trace  gases  are  separated  into 
the  radical  species  and  the  long-lived  species.  For  the 
radical  species,  local  photochemcial  equilibrium  is  assumed. 
Thus  the  concentrations  are  given  by  the  system  of  algebraic 
equations : 

Q.(N.)  -  0  i  «  1,  ...NR  where  NR  is  the  (A. 45) 

1  number  of  radical 

species 

j  *  1,  ...NS  where  NS  is  the 

total  number  of  trace 
species  in  the  model 


The  radical  species  treated  in  the  model  is  listed  in 
Table  A- 2  and  grouped  into  families. 

In  the  grouping  technique,  the  solution  to  Equation 
(A. 45)  provides  the  partition  of  the  radical  species  within 
each  family.  The  abundance  of  each  family  group  is  related 
to  long-lived  species.  The  abundance  of  the  oxygen  family 
is  to  a  large  extent  controlled  by  C>2  concentration,  while 
the  hydrogen  and  methyl  family  depends  on  water  vapor  con¬ 
centrations.  The  treatment  of  O2  and  H2O  concentration  is 
summarized  in  Table  A-3. 

The  sums  of  the  nitrogen  radicals  (NOX)  and  the 
chlorine  radicals  (C1X)  are  each  treated  as  a  long-lived 
species.  The  ClX  concentration  is  obtained  from  the  solution 
of  the  trace  gas  equation  with  photo-oxidation  of  CC14 , 

CFC13,  CF2C12,  CH^Cl  and  CH3CC13  as  the  source  and  tropo¬ 
spheric  washout  as  sink.  In  the  case  of  NOX,  photo-oxidation 
of  N20  constitutes  the  source  while  the  sink  consists  of 
tropospheric  washout  and  stratospheric  removal  via  the 
reaction. 


N  +  NO-*-  N2  +  0. 


The  production/loss  rates  of  the  long-lived  species 
are  expressed  in  the  form 


Qi 


pi  -  Liwfi  -  l1(2)  mL)2 


(2) 

where  is  the  production  rate,  and  are  the  linear 

and  quadratic  loss  frequency  respectively. 
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Table  A-2 

Group  of  Radical  Species 


Family 

Species 

Oxygen 

0,  0 ( ' D) 

Hydrogen 

H ,  OH,  H02,  H202 

Methyl  family 

ch3,  ch3o2,  ch3o,  ch2o,  hco 

Nitrogen 

NO,  N02,  N03,  N205 

hno3,  cino3,  n,  h^o4 

Chlorine 

Cl,  CIO,  C1N03,  HC1,  HOC1 
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Table  A- 3 


Fixed  Species 


_ Remark _ 

Calculated_from  ideal  gas 
law  using  T 

Set  equal  to  21%  of  N 


Stratospheric  value  of 
5  ppmv  is  assumed 

Tropospheric  mixing  ratio 
is  calculated  from  lati¬ 
tudinal  dependant  relative 
humidity  profile  (see  Table 
A-5). 


A  zero  flux  boundary  condition  is  imposed  at  the 
poles  and  the  top  t-undary  for  all  long-lived  species  except 
03  for  which  photochemical  equilibrium  is  assumed  at  the  top 
of  the  atmosphere.  Lower  boundary  conditions  are  specified 
to  simulate  appropriate  exchange  with  the  surface.  The  list 
of  lor "-lived  species  in  the  model  is  given  in  Table  A-4 
together  with  the  lower  boundary  conditions. 

The  reaction  rate  constants  are  taken  from  NASA/WMO 
workshop.  Tropospheric  washout  is  parameterized  by  a  first 
order  loss  frequency  with  effective  lifetime  of  10  days. 

A  24-hour  averaging  is  used  to  calculate  the  photolysis 
rate  to  simulate  diurnal  averaged  condition. 

A. 8  Model  Algorithm 

The  model  algorithm  employed  in  the  present  model  is 
schematically  represented  in  Figure  A-8.  For  our  numerical 
studies,  the  radical  species  concentrations  and  the  produc¬ 
tion  rates  and  loss  frequencies  of  the  long-lived  species  are 
calculated  once  every  10  days.  The  long-lived  species 
solutions  are  propagated  forward  in  time  for  10  days  with 
the  production  rates  and  loss  frequencies  held  fixed.  The 
values  for  the  eddy  diffusion  tensor  is  updated  once  a  month. 
The  temperature  and  1 1>  filed  is  read  in  and  interpolation 
for  each  time  step. 
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Table  A- 4 


Long-Lived  Species  and  Lower  Boundary  Condition 


i 

i. 


re 


►  ■  • 


e 


Species 

Lower  Boundary  Condition 

300  ppbv 

NOX 

. 1  ppbv 

ch3ci 

. 8  ppbv 

cci4 

.12  ppbv 

CH3CC13 

.12  ppbv 

CFC-11  7 

Time-dependent  release  rate 

CFC-12  r 

from  CMA  (1980)  and  spread 
uniformly  between  24 °N  and 
60°N 

C1X 

1  ppbv 

H2 

500  ppbv 

CO 

90°S-50°S  60  ppbv 

40°S-Eq.  70  ppbv 

Eq-40°N  140  ppbv 

80°N-90°N  190  ppbv 

CH4 

1.5  ppmv 

°3 

30  ppbv 
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Table  A- 5 


Relative  Humidity 


\Lat. 

Alt.  \ 

90S-45S 

42S-24S 

24S 

+ 

5S 

5S 

+ 

5M 

5N 

+ 

24N 

24N-42N 

42N-90N 

£  «  0 

.65 

.5 

.6 

.5 

.65 

i  *  0.5 

.5 

.38 

.43 

.38 

.5 

mss 

.55 

.58 

.55 

H§H 

§§§ 

.8 

m 

§||||| 

Stratospheric  value  of  5  ppmv  is  used. 


N 
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Figure  A- 8 

Schematic  Flow  Chart  for  2 


@ 


0  Model 
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APPENDIX  B 


LATITUDE-ALTITUDE  CROSS-SECTIONS 
FOR  RADICAL  SPECIES 


"1 


1 


N 


Presented  here  are  the  calculated  latitude-altitude 
cross-sections  for  some  of  the  radical  species  from  the  AER 
2-D  model  corresponding  to  March  conditions.  The  concentra¬ 
tions,  which  are  computed  using  24-hour  averaged  photolysis 
rates,  can  be  interpreted  as  the  diurnal  average  species 
concentrations.  This  should  be  taken  into  account  in  com¬ 
paring  the  calculated  results  with  available  measurements. 

The  behavior  of  the  species  in  the  upper  stratosphere 
can,  to  a  large  extent,  be  explained  in  terms  of  photochem¬ 
ical  interactions.  Some  of  the  species  exhibit  sizeable 
latitudinal  variations  in  the  lower  stratosphere.  In  the 
case  of  the  radical  species  with  longer  photochemical  life¬ 
times  (HC1 ,  HNO3) ,  the  effect  of  transport  is  evident. 
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